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I. 


INTRODUCTION 


For  over  150  years,  spectral  remote  sensing  has  been  an  acknowledged  utility, 
dating  back  to  the  invention  of  the  camera.  While  spectral  remote  sensing  has  enjoyed 
substantial  development  since  the  camera's  inception,  recent  advances  in  video,  scanners, 
cameras,  and  other  remote  sensing  support  systems  have  produced  sensor  platforms  of 
high  spatial,  spectral,  and  radiometric  resolution. 

In  September  1999,  the  Ikonos  satellite  was  launched  by  Lockheed  Martin  for 
Space  Imaging.  Ikonos  had  the  world's  highest  resolution  among  commercial  satellites 
for  two  years.  Its  camera,  built  by  Eastman  Kodak,  provides  ImlWbit  panchromatic 
and  Am  mbit  four  band  multispectral  data.  This  represents  a  substantial  improvement 
in  resolution  over  the  Landsat  7  satellite  launched  in  April  1999,  that  supports  the 
Enhanced  Thematic  Mapper  Plus  sensor. 

The  utility  of  image  classification  schemes  is  affected  by  increased  resolution.  In 
a  classified  image,  high  resolution  may  result  in  classified  regions  without  sufficiently 
defined  boundaries.  In  general,  in  some  instances,  increased  resolution  merely  introduces 
noise  that  obscures  the  features  of  most  interest  to  the  observer  (Wilkie  and  Einn,  1996). 
However,  the  use  of  appropriate  training  data  sets  in  performing  a  supervised  image 
classification  may  provide  the  necessary  means,  variances,  and  covariances  to  combat 
this  noise  associated  with  high  resolution. 

The  goal  of  this  thesis  is  to  create  a  class  map  of  the  Elkhom  Slough  and  its 
surrounding  wetlands.  The  Environment  for  Visualizing  Information  (ENVl)  software 
and  a  reference  class  map  of  the  study  area  provided  by  the  Elkhorn  Slough  Eoundation 
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will  be  used.  A  eomparison  of  the  created  and  reference  class  maps  will  be  made.  The 
comparison  will  be  used  to  evaluate  the  utility  of  using  four  band  high  resolution  data  and 
training  data  to  perform  scene  classification. 

Chapter  II  of  this  thesis  will  present  a  description  of  spectral  remote 
sensing  and  multispectral  imaging  theory.  Chapter  III  discusses  the  means  by  which 
multispectral  data  is  acquired  and  corrected.  Multispectral  data  and  classification 
schemes  are  explained  in  Chapter  IV.  The  Ikonos  and  Landsat  7  remote  sensing  systems 
are  described  in  Chapter  V.  Chapter  VI  presents  the  study  area  and  spectral  analysis. 
Chapter  VII  provides  a  summary  of  the  spectral  analysis  and  causes  of  any  inaccuracy. 
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II,  MULTISPECRTAL  IMAGING  THEORY 


A,  SPECTRAL  REMOTE  SENSING 

1,  Background 

Remote  sensing  is  the  proeess  of  aequiring  and  reeording  information  about  an 
objeet  without  being  in  direet  eontaet  with  that  objeet  (Gibson,  2000).  Although  this 
information  may  present  itself  in  various  forms,  sensing  the  electromagnetie  speetral 
information  produced  by  an  object  is  of  concern  here.  The  process  of  sensing 
electromagnetic  spectral  information  is  referred  to  as  spectral  remote  sensing.  The 
unique  way  that  each  object  emits,  reflects,  absorbs,  or  scatters  electromagnetic  radiation 
allows  an  electromagnetic  sensing  device  to  discriminate  between  different  objects. 
Objects  typically  produce  a  spectral  signature  by  reflecting  electromagnetic  energy 
produced  from  a  source  that  could  be  local  or  non-local  to  the  remote  sensing  platform 
that  contains  the  sensing  device.  These  two  processes  are  referred  to  as  active  and 
passive  remote  sensing  respectively.  The  sensor  collects  this  energy  in  spectral  bands 
defined  at  particular  regions  of  the  electromagnetic  spectrum.  The  energy  is  then 
converted  into  computer  processable  data  that  can  be  displayed  in  a  more  usable  form 
called  an  image. 

An  essential  element  in  remote  sensing  is  the  ability  to  record  images  of  the 
earth's  surface  (Campbell,  1996).  Space  borne  or  airborne  sensors  are  typically  used  to 
record  these  images  as  their  operational  altitudes  provide  a  means  for  collecting  spectral 
information  from  large  areas.  They  provide  us  with  as  much  or  more  detail  than  field 
surveys,  accomplish  this  over  large  areas,  and  can  do  so  repeatedly  (Wilkie  and  Finn, 
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1996).  The  spectral  characteristics  of  a  region  or  scene  can  provide  detailed  information 
about  that  scene  such  as  its  human  land  use,  water  resources,  and  distribution  of 
vegetation  as  multiple  images  can  be  obtained  simultaneously  from  different  wavelength 
bands.  Hence  using  multispectral  imaging,  spectral  information  from  the  near  infrared, 
short-wave  infrared,  mid-wave  infrared,  and  long-wave  infrared  region  can  be  collected. 
These  regions  are  outside  the  visible  portion  of  the  electromagnetic  spectrum.  These 
bands  are  used  to  discriminate  features  that  are  not  visible  to  the  human  eye 
(Multispectral  Users  Guide,  1995). 

Collecting  spectral  data  using  space  borne  or  airborne  sensors  is  a  nonintrusive 
way  to  obtain  information  that  could  provide  support  to  address  questions  concerning 
conserving  a  biological  system.  Human  subsistence  and  market  activities  are  bringing 
increasing  distress  to  ecological  systems  around  the  globe,  jeopardizing  the  structure, 
function,  production,  and  resilience  of  earth's  life  support  systems  (Wilkie  and  Finn, 
1996). 


2.  Electromagnetic  Spectrum 

The  electromagnetic  energy  collected  by  a  remote  sensing  device  can  be 
characterized  by  its  distribution  of  wavelengths.  A  continuum  consisting  of  all  possible 
wavelengths  is  called  the  electromagnetic  spectrum.  Figure  2.1  defines  some  useful 
regions  of  the  electromagnetic  spectrum.  However  the  range  of  wavelengths  between 
0.4 to  14.0^m  is  most  used  in  spectral  remote  sensing.  This  is  attributable  to  the 
existence  of  several  wavelength  sub-bands  between  0.4 and  14.0^m  that  allow 
partial  or  complete  transmission  of  electromagnetic  radiation  through  the  atmosphere. 
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These  wavelength  sub-bands  are  referred  to  as  atmospheric  windows.  The  atmosphere 
attenuates  radiation  passing  through  it  on  a  one-way  trip  for  about  100  kilometers.  If  a 
sensor  is  carried  by  a  low-flying  aircraft,  effects  of  the  atmosphere  on  image  quality  may 
be  negligible  (Campbell,  1996). 

Figure  2.2  shows  the  locations  of  atmospheric  windows  used  in  remote  sensing. 
Atmospheric  windows  are  defined  by  regions  of  high  transmission.  Passive  remote 
sensing  systems  take  advantage  of  this  by  using  the  sun's  illumination  to  sense  at 
wavelengths  from  0.4 to  3.0 //m.  Active  remote  sensing  systems  must  generate 
radiation  within  an  atmospheric  window.  The  most  common  active  remote  sensing 
system  is  radar  that  produces  electromagnetic  radiation  in  the  microwave  band  of  the 
spectrum  (Gibson,  2000). 

3.  Spectral  Response 

The  spectrum  of  reflected  wavelengths  for  reflective  materials  can  be  depicted 
graphically  as  a  reflectivity  curve.  The  reflectivity  curve  yields  reflectance  as  a 
percentage  of  incident  radiation  reflected  from  a  specific  material  as  a  function  of  the 
incident  radiation's  wavelength.  The  percentage  is  called  albedo.  Figure  2.3  presents 
typical  albedos  of  several  materials  observed  in  the  visible  and  near  infrared  bands. 

The  curves  differ  dramatically  as  each  material  interacts  with  the  sun's  radiation 
to  reflect  energy  in  distinct  ways.  For  healthy  vegetation,  chloroplasts  that  contain 
chlorophyll  reflect  wavelengths  corresponding  to  greenish  colors  (between  0.7  /ilm  and 
1.3  ^m)  as  the  chloroplasts  absorb  all  other  wavelengths.  In  the  process  of  absorption, 
the  incident  radiant  energy  is  transformed  and  not  reflected.  Between  10  and  30  percent 
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of  the  total  amount  of  visible  light  arriving  at  the  surface  of  a  leaf  is  reflected  back  as 
green  (Gibson,  2000).  The  near  infrared  wavelengths  yields  an  even  greater  spectral 
response  in  healthy  vegetation  as  mesophyll  cells  in  leaves  reflect  about  60  percent  of  the 
incident  near  infrared  radiation  that  reaches  the  leaves.  Since  the  abundance  of 
chlorophyll  and  mesophyll  in  plants  with  leaves  are  directly  linked  to  their  health,  poor 
plant  health  yields  lesser  reflectance  in  both  the  green  and  the  near  infrared. 

The  presence  of  rocks,  soils,  or  moisture  with  vegetation  is  almost  inevitable. 
This  will  significantly  affect  the  spectral  signature  of  a  scene  said  to  be  vegetation.  Soils 
constitute  the  primary  background  spectral  response  in  most  vegetated  landscapes,  and,  in 
combination  with  rock,  the  sole  spectral  response  in  nonvegetated  terrain  (Wilkie  and 
Finn,  1996).  Soils  from  different  regions  may  produce  different  spectral  responses  as  the 
mineral  composition  of  soils  varies  from  region  to  region.  Mineral  reflectance  spectra  are 
strongly  dependent  on  the  chemical  composition  of  the  rock  (Elachi,  1987).  Wavelengths 
from  2.08  /ilm  to  2.35  jilm  have  strong  absorptions  in  materials  containing  hydroxyl  ions 
making  these  wavelengths  useful  for  identifying  hydrothermally  altered  rocks. 

Moisture  tends  to  reduce  the  overall  reflectance  of  all  materials.  As  the  moisture 
content  of  a  material  increases,  the  material's  reflectance  tends  towards  the  reflectance  of 
water.  This  is  partly  responsible  for  enabling  one  to  distinguish  between  soil  regions 
differing  in  extents  of  irrigation.  Likewise,  one  is  also  able  to  distinguish  between  water 
regions  with  different  soil  concentrations.  For  example,  a  pond  is  a  body  of  still  water, 
and  under  most  conditions  will  have  a  different  soil  concentration  than  a  river,  stream,  or 
any  other  body  of  running  water.  In  these  cases  the  one  would  expect  the  reflectivities  to 
be  closer  to  clear  water  than  soil  or  even  moist  soil  but  still  discernable  from  one  another. 
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Clear  water  appears  blue  beeause  of  its  high  refleetanee  between  wavelengths  of  OA  /im 


and  0.6  /im  and  appears  blaek  at  wavelengths  beyond  0.7  /ilm .  Water  appears  blaek  in 
the  infrared  where  its  absorption  is  relatively  high  and  refleetanee  is  low.  Figure  2.4 
illustrates  how  deeply  radiation  ean  penetrate  into  clear  water.  When  soil  or  other 
particulate  is  suspended  in  water,  depending  on  the  depth  and  concentration  of  the  soil  or 
particulate,  the  spectral  reflectance  deviation  from  clear  water  can  be  dramatic.  For 
example,  plankton  reflects  in  the  green  and  near  infrared.  However,  water  containing 
plankton  only  reflects  in  the  green.  Figure  2.4  shows  that  the  plankton  filled  water  is 
green  because  green  wavelengths  not  only  penetrate  to  reach  more  plankton,  they  emerge 
from  the  water  with  far  less  attenuation  than  the  infrared  wavelengths.  This  idea 
becomes  important  when  viewing  water  in  a  false  color  composite. 

4.  General  Radiation  Laws 

To  be  modeled  as  a  black  body,  a  body  must  emit  the  maximum  intensity  of 
radiance  across  all  wavelengths  that  it  is  theoretically  possible  to  radiate  for  the  particular 
temperature  of  the  body  (Gibson,  2000).  True  black  bodies  do  not  exist;  a  black  body  is 
an  idealized  concept.  Being  a  black  body  requires  it  to  absorb  all  radiation  incident  upon 
it  and  emit  radiation  at  the  maximum  theoretical  rate.  In  defining  this  maximum 
theoretical  rate,  it  is  instructive  to  begin  with  the  production  of  radiation. 

a.  Electromagnetic  Radiation 

Electromagnetic  radiation  is  generated  by  the  conversion  of  energy  from 
other  forms  such  as  kinetic,  chemical,  or  nuclear  (Elachi,  1987).  Materials  generate  this 
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radiation  as  their  molecules  seek  lower  energy  states.  The  frequency  of  radiation 
generated  by  a  material  depends  upon  the  process  by  which  the  material's  molecules  seek 
lower  energy  states.  Electronic  transitions  in  atoms  tend  to  produce  radiation  in  the 
visible  and  ultraviolet  wavelengths.  Infrared  wavelengths  are  produced  by  vibrational 
and  rotational  transitions  in  molecules. 

When  radiation  is  produced  by  any  material,  it  may  possess  a  degree  of 
polarization.  Electromagnetic  radiation  consists  of  waves  with  alternating  electric  and 
magnetic  fields.  These  electric  and  magnetic  fields  are  perpendicular  to  each  other  and 
are  in  a  plane  perpendicular  to  the  direction  of  wave  propagation.  However  the 
orientation  of  the  fields  relative  to  a  coordinate  frame  may  be  defined  in  any  direction. 
The  orientation  of  the  fields  is  termed  the  polarization  of  the  radiation  (Rees,  1990).  The 
types  of  polarization  that  the  fields  could  have  are  called  plane,  circular,  elliptical,  and 
random  polarization.  In  plane  polarization  the  net  electric  and  magnetic  fields  each 
sweep  out  planes  as  the  wave  propagates.  In  circular  polarization  the  net  electric  and 
magnetic  each  fields  sweep  out  helices  as  the  wave  propagates,  and  the  helices  can  rotate 
either  clockwise  or  counter  clockwise  about  the  direction  of  propagation.  These  rotations 
further  define  circular  polarization  as  being  either  left-hand  circularly  polarized  or  right- 
hand  circularly  polarized.  A  simple  consistent  convention  would  distinguish  the  two.  In 
elliptical  polarization  the  net  electric  and  magnetic  fields  each  sweep  out  elliptical  helices 
as  the  wave  propagates.  In  random  polarization  the  radiation  is  a  mixture  of  the  plane, 
circularly,  and  elliptically  polarized  waves.  Radiation  in  this  condition  is  also  said  to  be 
unpolarized.  Sun  radiation  is  said  to  be  unpolarized  (Gibson,  2000).  The  importance  of 
polarization  lies  in  the  fact  that  some  sensing  devices  are  designed  to  only  detect 
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radiation  of  a  particular  polarization,  and  processes  of  refleetion  and  absorption  may  alter 
the  polarization  of  radiation.  The  effeets  of  polarization  generally  do  not  beeome 
considerable  until  sensors  try  to  colleet  radiation  at  relatively  long  wavelengths. 
Polarization  is  not  considered  in  the  visible  and  infrared  region,  but  is  eonsidered  in  the 
mierowave  region  used  by  radar. 

b.  Planck  Radiation  Laws 

It  was  mentioned  that  the  sun  is  a  source  of  polarized  radiation.  It  is  also  a 
source  that  eould  be  considered  a  blaek  body  as  previously  defined.  The  sun  has  a 
surfaee  temperature  of  5,750-6,000  K  and  radiates  across  a  range  of  wavelengths  at  an 
average  distance  from  the  Earth  of  150  million  kilometers  (Gibson,  2000).  To  compute 
the  distribution  of  energy  among  these  wavelengths  the  Planek  Radiation  Law  applies. 
This  radiation  law  is  Equation  2.1.  Equation  2.1  yields  the  speetral  emittanee  of  a  blaek 
body  as  a  eurve.  Eigure  2.5  in  Appendix  A  gives  the  speetral  emittanee  of  three  blaek 
body  sourees. 


^  - to -  [Equation  2.1] 

c^(e*^  - 1) 

h  =  Planek's  constant  ( 6 . 626 1  x  1 0“^"^  J -s) 

T  =  Temperature(  K  ) 

k  =  Boltzmann  constant  (1.3807xl0“^V- ) 
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c  =  Speed  of  light  (3x10*  —  ) 

s 

A,  =  Wavelength  (m) 

Like  the  sun,  the  earth  is  also  eonsidered  a  blaek  body.  The  peak  of  the 
eurve  represents  the  wavelength  that  eontains  the  greatest  intensity  of  emitted  radiation. 
To  compute  this  wavelength  the  Wien  Displacement  Law  applies.  This  law  is  Equation 

2.2. 


A  max  =  —  [Equation  2.2] 

5  =  2,S9S^m-K 

T  =  Temperature  {K) 

Eohowing  the  equation  development  of  W.G.Rees,  1990,  to  compute  the 
total  radiant  exitance  produced  by  a  black  body,  first  the  spectral  emittance  (Equation 
2.1)  is  integrated  over  ah  wavelengths  or  frequencies  to  yield  L.  It  is  then  used  in 
Equation  2.3  to  define  the  total  radiant  exitance  (Equation  2.4). 

Irradiance=  J  L(cosff)dQ.  [Equation  2.3] 

2k 

6  =  Angle  of  incidence  relative  to  the  normal 

JQ  =  Solid  angle  containing  the  incoming  radiation 
L  =  Total  radiance 
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[Equation  2.4] 


—  iTT^k 

M  =  AT'^  where  A  = - ^ 

\5c^h^ 


The  A  is  called  the  Stefan-Boltzmann  constant  and  its  value  is  5.61'K\0~^Wm~^K~^ . 
These  equations  become  very  useful  calibrating  remote  sensing  devices. 

Recall  that  Equation  2.1  and  Equation  2.4  are  based  upon  the  idea  that  the 
materials  whose  spectral  characteristics  they  describe  are  ideal  sources  and  sinks.  Real 
sources  and  sinks  involve  less  than  100%  radiative  processes.  These  processes  are 
absorption,  reflection,  and  scattering.  They  are  measured  as  ratios  of  their  individual 
radiant  powers  to  incident  radiant  power.  If  ^  is  the  absorptance,  p  is  the  reflectance, 
and  T  is  the  transmittance,  then  Equation  2.5  must  hold. 

P  +  p  +  T  =  \  [Equation  2.5] 


If  an  object  has  the  properties  (3  =  1  and  /7  =  T  =  0,  then  the  object  is  a  black  body 
(Seyrafi  and  Hovanessian,  1993).  Equations  2.1  and  2.4  must  be  modified  to  describe 
real  sources  and  sinks.  The  modification  factor  is  called  emissivity.  Kirchoff s  law  states 
that  the  ratio  of  emitted  radiation  flux  to  absorbed  radiation  flux  is  the  same  for  all 
blackbodies  at  the  same  temperature  (Campbell,  1996).  This  ratio  is  emissivity.  It  is 
defined  by  Equation  2.6. 


4. 


[Equation  2.6] 


II 


L^P  =  Black  Body  Radiance  at  wavelength  A 


p  =  l-£  [Equation  2.7] 
p  =  Reflectivity  with  zero  transmittance 

Materials  with  £  =  \  are  black  bodies,  with  f  =  0  are  called  white  bodies,  and  with 
0<£'<1  are  called  gray  bodies.  Equations  2.1  and  2.4  are  modifled  by  simply 
multiplying  them  by  the  appropriate  emissivity  to  describe  real  bodies. 

Emissivity  becomes  very  important  when  sensing  radiation  from  7  pm  to  1 8  pm , 
the  thermal  infrared  region.  Many  of  the  same  kinds  of  film,  filters,  and  cameras  that  we 
use  in  the  visible  portion  of  the  spectrum  can  also  be  used,  with  minor  variations,  for 
imaging  in  the  near  infrared  (Campbell,  1996). 

c.  Spatial  Resolution  Limitations 

Diffraction  is  another  process  involving  electromagnetic  radiation  of 
concern  in  spectral  remote  sensing  because  of  its  importance  in  establishing  the  optical 
resolving  power  of  remote  sensors.  W.G  Rees,  1990,  shows  that  the  amplitude  of  the 
diffracted  field  in  the  case  of  Eraunhofer  diffraction  (far-field  diffraction)  for  a  circular 
aperture  is  given  by  Equation  2.8. 


Amplitude  =  Q  [Equation  2.8] 

{KDsmOJl) 
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Jj  =  First  order  Bessel  function 

K  =  Wave  number 
D  =  Diameter  of  circular  aperture 
6^  =  Radial  angle 
Q  =  Proportionality  Constant 

Figure  2.6  illustrates  the  Fraunhofer  diffraction  geometry  for  a  circular  aperture.  The 
locations  of  the  zeros  are  computed  using  Equation  2.9. 

sin(^^)  =  1.22^  [Equation  2.9] 

A.  =  Wavelength 

A  useful  criterion  for  establishing  the  spatial  resolution  limit  is  that  the 
maxima  of  the  amplitude  produced  by  one  object  must  just  coincide  with  the  first 
vanishing  of  the  amplitude  produced  by  the  other  object.  The  angular  separation 
corresponding  this  diffraction  limit  or  ultimate  resolving  power  is  defined  by  Equation 

2.9. 

In  all  optical  systems  diffraction  sets  the  ultimate  resolution  limit, 
although  it  may  not  be  responsible  for  the  operational  resolution  limit.  Eimitations 
provided  by  other  system  factors  can  be  more  restrictive  on  resolution  than  diffraction. 
The  resolution  of  an  optical  system  is  a  combination  of  the  performances  of  the  lens  and 
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the  film  (Rees,  1990).  When  a  digital  system  is  used,  resolution  is  limited  by  the 
instantaneous  field  of  view  (IFOV). 

The  Modulation  Transfer  Function  (MTF)  is  used  to  measure 
performances  of  lenses.  The  MTF  for  a  given  spatial  frequency  expresses  the  fraction  of 
an  incoming  signal  of  that  frequency  that  the  detector  outputs  ("Figures  of  Merit"). 

B,  SENSING  THE  EARTH’S  SURFACE 

1,  Radiation  Propagation  in  the  Atmosphere 

Radiation  from  the  sun  has  to  propagate  through  the  atmosphere  twice  before 
being  collected  at  a  space  borne  sensor.  Solar  radiation  is  absorbed,  scattered,  and 
refracted  twice  before  reaching  the  sensor.  The  extent  of  the  atmosphere  is  about 
100A:mand  Table  2.1  gives  its  gaseous  composition.  COj  levels  have  become  of 
particular  concern.  COj  has  strong  absorption  in  the  spectral  region  from  lA/im 
\o\l  lim .  This  region  covers  a  large  portion  of  the  earth's  spectral  emittance.  Hence  COj 

absorbs  a  large  portion  of  the  radiation  emitted  by  the  earth. 

In  addition  to  the  gases  in  Table  2.1,  the  atmosphere  contains  water  vapor, 
methane,  and  aerosols  such  as  dust,  volcanic  ash,  smoke,  and  pollen.  The  two  main 
causes  of  attenuation  in  the  atmosphere  are  molecular  absorption  by  several  minor 
constituents  of  the  atmosphere  and  scattering  due  to  aerosols  (Seyrafi  and  Hovanessian, 
1993).  Figure  2.7  illustrates  the  spectral  irradiance  at  the  earth's  surface  and  outside  the 
earth's  atmosphere.  Take  note  of  the  regions  defining  atmospheric  windows.  Absorption 
is  primarily  due  to  electronic  transitions  in  atoms  and  vibrational  and  rotational 
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transitions  in  molecules.  As  mentioned  previously,  these  processes  produce  strong 
absorption  in  the  visible,  ultraviolet,  and  infrared  regions. 

Considering  only  the  gases  listed  in  Table  2.1,  attenuation  is  dominated  by 
Rayleigh  scattering.  Following  Seyrafi  and  Hovanessian,  1993,  Rayleigh  scattering 
attenuation  is  computed  using  Equation  2.10. 

/  =  [Equation  2.10] 

=  Original  Radiation 
/  =  Scattered  Radiation 
<7^  =  Rayleigh  scattering  coefficient 
R  =  Path  length 


The  Rayleigh  cross  section  presented  to  an  incident  wave  is  defined  so  that  the 

total  power  scattered  by  a  gaseous  molecule  in  all  directions  is  equal  to  the  power 
flowing  across  the  cross  section  (Seyrafi  and  Hovanessian,  1993).  This  cross  section  is 
computed  using  Equation  2.11. 


8;r'(n'-l) 

3AU' 


[Equation  2.1 1] 


n  =  Index  of  refraction  of  the  gaseous  medium 
N  =  Molecular  number  density 


15 


When  considering  only  the  aerosols,  attenuation  is  dominated  by  Mie  Scattering. 
Again  following  Seyrafi  and  Hovanessian,  1993,  Mie  scattering  attenuation  is  computed 
using  Equation  2.12. 


I  =  I„e  where  <J  (h)  =  <7  (0)  [Equation  2.12] 

M(0) 

=  Intensity  of  original  Radiation 
/  =  Intensity  of  scattered  Radiation 
(j^  =  Rayleigh  scattering  coefficient 
R  =  Path  length 

M(/z)  =  Aerosol  number  density  at  height  h 
<7^  (0)  =  Aerosol  coefficient  at  sea  level 

Eigure  2.8  gives  a  representative  profile  of  the  aerosol  density  versus  altitude. 

An  accurate  description  of  the  combined  effects  atmospheric  attenuations  on 
radiation  propagation  is  sufficiently  complex  to  require  computer  models.  Some  widely 
used  models  include  HITRAN  and  MODTRAN. 

Refraction  alters  the  direction  of  propagation  of  radiation  as  is  travels  from  one 
medium  to  another  with  a  different  index  of  refraction.  Snell's  Eaw  (Equation  2.13) 
involves  the  indices  of  refraction  of  both  mediums  describing  this  change  in  the  direction 
of  propagation.  The  index  of  refraction  is  defined  as  the  ratio  between  the  velocity  of 
light  in  a  vacuum  (c)  to  its  velocity  in  the  medium  (c„)  (Campbell,  1996).  Equation 
2.14  defines  the  index  of  refraction. 
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nsin(^)  =  n  sin(0  )  [Equation  2.13] 

=  Index  of  refraction  and  incident  angle  of  first  medium 

I  I 

n  ,0  =  Index  of  refraction  and  transmission  angle  of  second 
medium 

n  =  —  [Equation  2.14] 

Cn 

c„  =  Speed  of  light  in  material  with  index  of  refraction  n 

The  index  of  refraction  is  a  property  of  a  medium  and  is  determined  by  the  way  in 
which  atoms  respond  to  an  electromagnetic  wave.  This  response  can  be  influenced  by 
changes  in  the  thermodynamic  conditions  imposed  upon  the  medium,  such  as  temperature 
and  density.  The  atmosphere  has  layers  of  different  compositions,  temperatures,  and 
densities  making  refraction  a  considerable  effect  in  spectral  remote  sensing. 

Atmospheric  turbulence  is  characterized  by  random  air  motion.  This  motion 
causes  random  changes  in  the  local  index  of  refraction  bending  the  radiation  rays 
propagating  through  the  atmosphere.  It  is  the  cause  of  phenomena  such  as  the  twinkling 
of  stars  at  night,  shimmering  of  distance  objects  on  a  hot  day,  and  mirages  (Seyrafi  and 
Hovanessian,  1993).  Temperature  fluctuations  have  the  greatest  influence  on  the  index  of 
refraction  fluctuations.  Eigure  2.9  very  generally  illustrates  how  large  turbulent  wind 
eddies  can  form  to  cause  a  radiation  ray  to  wander. 
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The  effects  of  atmospheric  turbulence  can  be  characterized  as  changes  in  the 
beam  direction  (beam  wander),  changes  in  the  beam  size  (beam  breathing),  and  intensity 
fluctuations  (scintillation).  These  effects  manifest  themselves  as  motion,  blur,  and 
scintillation  in  the  image.  The  single  most  important  parameter  that  characterizes  the 
effects  of  turbulence  on  radiative  transfer  is  the  atmospheric  refractive  index  structure 
parameter,  C\. 

cl  is  computed  using  Equation  2.15. 


d  =  [79x10"'  [Equation  2.15] 

cl  =  Temperature  Structure  Parameter 

P  =  Atmospheric  pressure  ( mbar  ) 

T  =  Temperature  {K) 

Equation  2.16  is  given  by  Hufnagel  and  Stanley,  1964,  as  an  altitude-only  dependent  way 
to  compute  Cl . 


_i 

Cl  =  4.2xl0"'"/z'3e  [Equation  2.16] 

H  =  Altitude  above  the  ground  (m) 

=  3200m 
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Atmospheric  turbulence  can  reduce  the  performance  of  optical  imaging  systems, 
particularly  the  system's  resolution.  A  measure  of  the  extent  to  which  an  optical  system's 
performance  has  been  reduced  is  called  the  optical  transfer  function.  A  model  created  by 
D.L  Fried,  1966,  for  spherical  waves  is  given  by  Equation  2.17.  This  equation  is  valid 
only  for  exposures  taken  over  either  a  long  (long  term  OTF)  or  short  (short  term  OTF) 
time  interval. 


-i  -  fji  - 
-{21.62C^RA  V/ll-A—y]} 

OTF  =  e  [Equation  2.17] 

R  =  Path  length 

/I  =  Wavelength  of  the  radiation 
=  Angular  spatial  frequency 
=  Diameter  of  the  aperture 

A=  0  (Fong  term  OTF),  ^  (far-field  short  term  OTF),  1  (near¬ 
field  short  term  OTF) 


2.  Radiation  Interacting  With  Surface  Materials 

Radiation  interacts  with  surface  materials  through  the  processes  of  reflection, 
absorption,  scattering,  transmission,  and  fluorescence.  If  a  surface  is  smooth,  specular 
reflection  occurs  in  which  all  or  almost  all  incident  radiation  is  reflected  at  the  angle  of 
incidence.  If  the  surface  is  rough,  diffuse  reflection  occurs  in  which  all  or  almost  all 
incident  radiation  is  scattered  in  all  possible  directions  of  reflection.  A  perfectly  diffuse 
reflector  is  often  designated  as  a  Fambertian  surface  (Campbell,  1996).  Some  materials 


19 


have  a  complex  distribution  of  surfaces  that  receive  incident  radiation.  These  materials 
may  produce  specular  reflection  for  radiation  incident  at  one  angle  and  diffuse  reflection 
when  radiation  is  incident  at  a  different  angle.  A  bi-directional  reflectance  distribution 
function  (BRDF)  is  used  to  describe  the  complex  reflection  produced  by  a  surface  over 
all  possible  angles  incidence  and  reflection. 

Transmission  occurs  when  incident  radiation  is  transmitted  through  a  material 
with  little  or  no  attenuation.  The  transmittance  in  the  form  of  depth  penetration  of  water 
is  illustrated  in  Figure  2.4.  Even  though  a  material  may  reflect  radiation  of  a  particular 
wavelength,  penetration  depth  may  limit  the  extent  to  which  it  can  reflect,  thereby 
allowing  more  radiation  to  be  transmitted.  An  example  of  this  is  plant  leaves.  Plant 
leaves  reflect  infrared  radiation,  but  because  of  their  thin  structure  the  radiation  does  not 
penetrate  long  enough  to  be  completely  absorbed,  resulting  in  significant  amounts  of 
infrared  radiation  passing  through. 

Fluorescence  occurs  when  an  object  illuminated  with  radiation  of  one  wavelength 
emits  radiation  of  a  different  wavelength  (Campbell,  1996).  When  examining  remotely 
sensed  data,  one  needs  to  be  aware  that  there  are  some  types  of  minerals  and  plants  that 
fluoresce. 

The  combined  effects  of  reflection,  absorption,  and  transmission  are  defined  in 
Equation  2.18. 


Ej  (A)  =  (A)  -I-  E^  (A)  -I-  Ej.  (A)  [Equation  2.18] 

Ej  (A)  =  Incident  Energy 
(A)  =  Reflected  Energy 
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(A)  =  Absorbed  Energy 
Ej.  (A)  =  Transmitted  Energy 

Equation  2.18  is  simply  the  conservation  of  energy.  Reflectance,  transmittance, 
and  absorptance  are  as  previously  defined. 
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III.  COLLECTING  SPECTRAL  INFORMATION 


A,  Remote  Sensing  Platforms 

1,  Altitudes  of  Operation 

As  mentioned  previously,  speetral  data  is  typieally  eolleeted  utilizing  airborne  or 
spaee  borne  sensors.  This  is  beeause  speetral  data  ean  be  eolleeted  from  large  seenes  at 
high  altitudes.  Airborne  sensors  are  earned  by  aireraft  and  spaee  borne  sensors  are 
earned  by  satellites,  all  of  whieh  operate  at  various  altitudes.  Operational  altitudes  in  part 
determine  spatial  resolution  and  seale,  and  henee  data  utility.  Varying  the  distanee  at 
whieh  images  are  aequired  allows  different  types  of  information  to  be  extraeted  (Gibson, 
2000).  Figure  3.1  in  Appendix  A  illustrates  the  various  heights  at  whieh  typieal  remote 
sensing  deviees  operate. 

2,  Airborne  Sensors 

Airborne  sensors  basieally  eonsist  of  a  eamera  mounted  at  some  point  on  an 
aireraft.  The  eamera  may  obtain  images  at  an  oblique  angle  to  the  vertieal  or  at  the 
vertieal.  These  images  typieally  have  very  good  spatial  resolution  but  extensive 
geometrie  distortions.  The  geometrieal  distortions  often  manifest  themselves  as  ehanges 
in  seale  from  one  region  in  the  image  to  another.  Images  taken  at  the  vertieal  may  be 
used  to  produee  maps  only  after  the  images  have  been  proeessed  to  get  rid  of  geometrieal 
errors. 
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Three  types  of  cameras  often  used  as  airborne  sensors  are  metric,  strip,  and 
panoramic  cameras.  Although  these  cameras  vary  in  design,  some  of  the  basic 
components  are  the  same. 

a.  Camera  Components 

In  general,  light  enters  the  camera  through  the  lens  cone  assembly,  which 
contains  the  lens,  filter,  shutter,  and  diaphragm.  Passing  through  the  lens  cone,  light  is 
fdtered  to  the  desired  wavelength,  focused,  and  cut  for  exposure  time  and  space  on  film. 
A  film  magazine  holds  the  film.  It  usually  includes  a  supply  spool,  holding  perhaps 
several  hundred  feet  of  unexposed  aerial  film,  and  a  take-up  spool  to  accept  exposed  fdm 
(Campbell,  1996).  An  electrical  film  drive  device  progresses  the  film.  Cameras  use 
filters  to  collect  single  wavelength  band  imagery. 

b.  Strip  Cameras 

In  strip  cameras  the  film  is  passed  over  a  slit.  The  slit  is  analogous  to  a 
shutter  that  never  closes  from  the  start  to  the  end  of  image  collection.  Images  are  not 
separated  by  frames.  Once  image  collection  is  complete,  the  result  is  one  long  strip  of 
imagery.  Figure  3.2  illustrates  this  idea. 

The  speed  of  film  movement  as  it  passes  the  slit  is  coordinated  with  the 
speed  and  altitude  of  the  aircraft  to  provide  proper  exposure  (Campbell,  1996).  Strip 
cameras  provide  good  use  for  measuring  long  straight  regions.  Disadvantages  include 
image  distortion  due  to  any  change  in  the  aircraft's  motion. 
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c. 


Panoramic  Cameras 


Panoramic  cameras  use  a  rotating  lens  or  prism  to  direet  light  to  the  film. 
It  also  uses  a  slit  for  film  exposure,  however  the  slit  moves  with  a  prism  or  lens  aeross  the 
film  while  the  film  is  stationary.  The  film  remains  stationary  during  eaeh  sweep  of  the 
prism  and  lens.  This  proeess  results  in  images  consisting  strips  covering  wide  areas. 
Figure  3.3  illustrates  this  idea. 

Disadvantages  inelude  image  distortion  due  to  forward  motion  of  the 
aireraft  and  the  side-to-side  sean  of  the  lens  and  prism.  Panoramie  photographs  have 
serious  geometrie  distortions  that  require  eorreetion  before  they  ean  be  used  as  the  basis 
for  measurements  (Campbell,  1996).  The  prineiple  advantages  of  the  panoramie  eamera 
are  its  high  image  resolution  and  large  area  of  eoverage  (Lihesand  and  Kiefer,  1987). 

d.  Metric  Cameras 

Metrie  eameras  provide  images  with  the  least  geometrie  errors.  They  are 
often  referred  to  as  mapping  eameras.  A  low  distortion  lens  system  is  held  in  a  fixed 
position  relative  to  the  plane  of  the  film  (Lihesand  Kiefer,  1987).  Eaeh  time  the  shutter  is 
opened,  only  one  frame  of  imagery  is  obtained.  An  eleetronie  system  regulates  how  often 
the  shutter  opens. 


e.  Image  Distortions 

Images  are  primarily  taken  with  the  eamera  oriented  with  a  vertieal  view 
of  the  ground.  Three  types  of  geometrieal  errors  ean  be  charaeterized  as  terrain,  optieal, 
and  radial  displaeement  errors.  Optieal  distortions  are  eaused  by  defeets  in  the  lens. 
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Radial  displacement  oeeurs  when  objeets  appear  to  lean  outward  from  the  eentral 
perspeetive  of  the  eamera  lens  (Campbell,  1996).  Terrain  errors  are  produeed  as  ehanges 
in  the  elevation  of  a  seene  ehanges  the  seale. 

A  refining  of  lens  making  proeesses  has  made  optieal  distortions 
insignifieant.  Terrain  and  radial  errors  are  correeted  in  orthoreetifieation  or  rubber 
sheeting.  Rubber  sheeting  uses  a  mathematieal  model  to  adjust  the  image  aeeording  to 
ground  eontrol  points  (GCP).  Orthoreetifioation  is  a  proeedure  whieh  models  the 
viewing  geometry  of  the  sensor  and  aeeount  for  distortions  introdueed  by  terrain  relief 
("Remote  Sensing:  Key  Coneepts").  This  proeess  involves  seanning  the  image  into  a 
eomputer  system,  thereby  digitizing  it.  A  eomputer  program  then  uses  the  image,  with  a 
map  that  illustrates  the  eorreet  loeations  of  figures  in  the  image,  to  eorrect  the  image.  In 
order  for  this  transformation  to  be  aeeurate,  it  is  also  neeessary  to  input  a  digital  elevation 
model,  whieh  stores  the  height  of  all  the  seanned  pixels  (Gibson,  2000).  The  result 
orthoreetifieation  is  an  image  with  seale  variations  removed.  The  x  and  y  axes  of  the 
image  are  aligned  with  the  east  and  north  axes  of  the  map  projeetion  being  used,  and  eaeh 
pixel  in  the  image  is  geo-coded  with  its  assoeiated  coordinate  in  the  map  projection 
("Orthoreetifieation:  What,  Why,  How,  Who?"). 

f  Geographic  Information  Systems 

It  is  often  eonvenient  to  use  overlays  to  define  or  interpret  an  image.  An 
overlay  eonsists  of  image  data  of  partieular  geographie  qualities;  for  example  vegetation 
elasses  or  human  development  elasses.  Many  overlays  defining  earth  resourees  in  a 
seene  is  ealled  a  Geographie  Information  System.  Manual  geographie  information 
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systems  usually  comprise  several  data  elements  including  maps,  sheets  of  transparent 
materials  used  as  overlays,  aerial  and  ground  photographs,  statistical  reports  and  field 
survey  reports  (Star  and  Estes,  1990).  Overlays  are  able  to  be  integrated  into  a  single 
data  set  that  can  be  analyzed  using  a  computer. 

2,  Space  borne  Multispectral  Scanning  Systems 

Space  borne  Multispectral  scanning  systems  (MSS)  collect  spectral  information  in 
two  distinct  ways.  The  first  way,  which  is  called  traverse  scanning,  MSS  use  a  rotating 
mirror  to  record  a  swath  of  ground  perpendicular  to  the  satellite’s  path.  A  telescope 
collects  light  and  focuses  it  onto  fiber-optic  bundles  that  transmit  the  light  to  detectors. 
The  light  is  then  divided  into  a  number  of  wavelength  bands  corresponding  to  the 
system's  spectral  resolution  and  sent  to  the  detectors. 

The  instantaneous  field  of  view  (IFOV)  of  a  scanning  instrument  can  be  defined 
as  the  ground  area  viewed  by  the  sensor  at  a  given  instant  in  time  (Campbell,  1996)  The 
spatial  resolution  of  the  system  sets  the  IFOV. 

The  second  way  in  which  the  MSS  collects  spectral  data  is  employing  pushbroom 
scanning  and  using  a  Charged  Coupled  Imager  (CCI)  as  a  sensor.  This  system  does  not 
use  a  scanning  mirror,  but  relies  on  quantum  mechanical  shifting  that  occurs  within  the 
CCI. 


a.  Charge  Coupled  Imager 

A  CCI  consists  of  an  array  of  Charge  Coupled  Devices  (CCD's).  Fight 
entering  the  MSS  is  formed  into  a  beam  using  lenses  and  is  split  into  a  number  of  spectral 
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regions  with  diffraction  gratings  or  filters.  Each  CCD  is  made  out  of  metal-oxide 
semiconductor  materials  and  develops  a  quantum  mechanical  potential  well.  When  light 
strikes  the  CCD  a  number  of  charges  are  generated  proportional  to  the  intensity  of  the 
light.  The  collected  charges  are  shifted  down  the  array  and  are  converted  to  equivalent 
current  or  voltage  at  the  output  terminal  (Karim,  1900).  The  shifting  is  controlled  by  a 
clock,  allowing  scanning  to  be  done  electronically.  The  CCD's  are  sensitive  to  the  visible 
and  near  infrared  regions.  Responsivity,  detectivity,  and  quantum  efficiency  are  among 
various  parameters  used  to  quantify  the  performance  of  semiconductor  sensors. 

b.  Traverse  and  Pushbroom  Scanning 

Figures  3.4  and  3.5  demonstrate  the  pushbroom  and  traverse  scanning 

methods. 


3.  Other  Remote  Sensing  Systems 

In  addition  to  the  MSS,  there  are  the  Hyperspectral  Scanning  System  (HSS)  and 
the  Ultraspectral  Scanning  System  (USS).  In  these  systems,  spectral  data  is  collected 
simultaneously  from  many  narrower  wavelength  bands  than  in  the  MSS.  Figure  3.4 
illustrates  the  relative  difference  in  the  number  and  width  of  the  wavelength  bands  in 
which  each  system  operates. 

B,  RESOLUTION 

There  are  four  types  of  resolution  of  interest  in  remote  sensing.  They  are  spectral, 
temporal,  radiometric,  and  spatial  resolution.  Spectral  resolution  is  a  measure  of  how 
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well  a  sensor  can  distinguish  a  material  based  on  the  sensor's  number  and  widths  of 
wavelength  bands.  High  spectral  resolution  would  be  characteristic  of  a  system  using  a 
large  number  of  small  wavelength  bands.  Temporal  resolution  is  a  measure  of  how  often 
data  can  be  obtained  from  the  same  area  (Gibson,  2000).  Spatial  resolution  is  a  measure 
of  the  smallest  spatial  length  a  material  can  possess  to  be  resolved.  Refer  to  Chapter  II 
for  detail  regarding  spatial  resolution.  Radiometric  resolution  measures  the  number  of 
shades  of  gray  a  sensor  can  detect  and  applies  only  to  digital  systems.  It  is  often 
expressed  in  bits  .  Equation  3.1  expresses  this  idea. 

Number  of  gray  levels  =  [Equation  3.1] 
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IV.  MULTISPECTRAL  DATA  ANALYSIS  TOOLS 


A,  MULTISPECTRAL  DATA 

Multispectral  data  can  be  imagined  to  be  a  number  of  registered  two  dimensional 
arrays  of  digital  numbers.  Eaeh  array  eorresponds  to  data  in  eaeh  speetral  band,  and  the 
digital  numbers  are  brightness  values  generated  by  the  remote  sensing  system.  Figure  4.1 
illustrates  this  for  remote  sensing  systems  of  different  speetral  resolutions. 

Figure  4.2  illustrates  the  digital  numbers  that  eorrespond  to  a  pixel.  Fleetronie 
sensors  reeord  digital  numbers  as  bits  in  the  binary  number  system.  The  number  of 
brightness  within  a  digital  image  is  determined  by  the  number  of  bits  available 
(Campbell,  1996). 

1.  Digital  Images 

In  a  given  band,  the  radianees  of  all  materials  eovered  by  a  pixel  eontribute  to 
determining  the  pixel's  digital  number.  A  digital  image  is  eomposed  of  an  array  of  pixels, 
with  eaeh  having  one  digital  number  per  band. 

When  examining  multispeetral  data  on  a  video  monitor,  the  digital  data  undergoes 
a  eomplex  proeess  in  whieh  band  brightness  is  eorresponded  with  portions  of  the  red, 
green,  or  blue  projeeted  on  the  sereen  by  the  monitor's  eleetronies.  It  is  often  useful  to 
examine  images  in  true  eolor  or  in  false  eolor.  Using  true  eolor,  the  red,  green,  and  blue 
bands  of  multispeetral  data  are  projeeted  onto  the  video  monitor  in  red,  green,  and  blue 
respeetively.  Using  false  eolor,  the  red,  infrared,  and  green  bands  are  projeeted  in  green, 
red,  and  blue  respeetively.  An  analyst  of  multispeetral  data  determines  to  use  either  a 
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true  color  or  false  color  image  based  on  how  well  the  image  distinguishes  important 
features  of  a  scene  under  investigation. 


B,  PRINCIPLE  COMPONENTS  ANALYSIS 

It  is  often  very  useful  to  plot  pixels  in  digital  number  space.  Digital  number  space 
is  a  multi-dimensional  space  in  which  each  axis  corresponds  to  each  available  band. 
Multispectral  imagery  tends  to  show  high  correlation  between  certain  bands  yielding 
information  redundancy.  A  principle  components  analysis  transformation  (PCA)  is 
performed  on  the  data  to  remove  this  correlation.  The  PCA  uses  a  variance-covariance 
matrix  to  calculate  a  coefficient  matrix,  which  is  applied  to  the  brightness  values  of  raw 
image  pixels  giving  a  decorrelation  (Wilkie  and  Finn,  1996). 

The  PCA  transform  is  a  linear  transformation  that  is  also  known  as  the  Karhunen- 
Loeve  or  Hotelling  transform  (Richards,  1995).  Equation  4.1  illustrates  the  PCA 
transformation.  Equation  4.1  gives  the  pixel  sets  before  and  after  a  PCA  transformation 
Eigures  4.3a  and  4.3b  illustrates  that  the  PCA  transformation  is  simply  a  rotation  in 
digital  number  space. 


Yi  =  GXi  [Equation  4.1] 

Yi  =  Vector  of  principle  component  values  of  the  pixel 

G  =  Einear  PCA  transform 

Xi  =  Vector  of  pixel  points  in  original  space 
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Figure  4.3b  illustrates  that  the  PCA  transformation  results  in  a  number  of  PCA 
eomponents.  Eaeh  eomponent  ean  be  represented  as  an  image.  The  number  of 
eomponents  produeed  by  the  transformation  equals  the  number  of  bands  in  the  raw  data. 
The  eomponents  posses  a  progressively  smaller  varianee.  The  first  three  eomponents  of 
Landsat  TM  imagery  will  typieally  eontain  over  95%  of  the  varianee  in  all  speetral  bands 
(Wilkie  and  Finn,  1996).  PCA  eomponents  beyond  the  first  few  typieally  show  the  noise 
from  the  original  data. 


C.  MULTISPECTRAL  IMAGE  CLASSIFICATION 

Informational  elasses  are  Regions  of  Interest  to  multispeetral  data  users.  These 
regions  may  be  forests,  areas  of  human  development,  or  any  other  region  to  be  analyzed. 
Speetral  elasses  are  groups  of  pixels  that  have  similar  speetral  properties.  When 
visualized  in  digital  number  spaee,  speetral  elasses  tend  to  form  elusters.  It  is  not 
uneommon  to  find  several  speetral  elasses  within  one  informational  elass.  For  example, 
an  informational  elass  defined  as  forest  may  possess  speetral  elasses  of  trees,  soil,  or  even 
shaded  areas.  The  goal  of  digital  image  elassifieation  is  to  assign  pixels  to  elasses  in  a 
region  under  investigation.  Eaeh  pixel  is  eompared  to  all  others  and  those  that  belong  to 
informational  elasses.  Classifieation  is  performed  using  a  eomputer  that  implements  any 
one  of  a  number  of  user  speeified  algorithms.  The  two  general  elassifieation  sehemes  are 
Supervised  and  Unsupervised  elassifieation.  This  thesis  will  only  present  some  detail 
regarding  supervised  elassifieation;  and  detail  regarding  unsupervised  elassifieation  may 
be  found  from  (Rees,  1996)  and  others. 
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1.  Supervised  Classification 

Supervised  elassifieation  is  the  proeess  of  assigning  unidentified  pixels  to 
informational  elasses  using  training  areas.  Training  areas  are  user  defined  samples  of 
pixels  with  known  identities.  Although  there  are  a  number  of  speeifie  algorithms  for 
performing  a  supervised  elassifieation,  this  thesis  will  offer  only  a  detailed  deseription  of 
the  Speetral  Angle  Mapper  and  the  Maximum  Likelihood  algorithms.  For  detail  of  other 
supervised  elassifieation  algorithms  see  Campbell,  (1996)  and  others. 

a.  Maximum  Likelihood 

The  maximum  likelihood  algorithm  is  eonsidered  to  be  the  most  eommon 
supervised  elassifieation  method  (Riehards,  1999).  This  elassifier  eonsiders  the  relative 
likelihood  that  a  pixel  belongs  to  a  partieular  class.  Probabilities  are  computed  from  the 
estimated  means  and  variances  of  the  classes.  A  pixel  is  then  assigned  to  the  class  that 
maximizes  the  probability.  Since  remotely  sensed  images  contain  classes  that  are  most 
often  impure,  the  maximum  likelihood  classifier  provides  the  advantage  of  taking 
variations  in  brightness  into  consideration.  Another  advantage  is  that  pixels  can  be 
classified  even  if  they  lie  within  an  overlap  region  defined  where  frequency  distributions 
of  different  classes  overlap.  Computation  of  the  estimated  probabilities  is  based  on  the 
assumption  that  both  training  data  and  the  classes  themselves  display  multivariate 
Gaussian  frequency  distributions  (Campbell,  1996). 

Before  the  classifier  proceeds  with  the  classification,  the  user  may  set 
probability  thresholds.  All  pixels  with  computed  probabilities  below  this  threshold  are 
unclassified.  The  results  of  this  classifier  are  decision  rules  and  a  classified  image.  A 
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decision  rule  is  created  for  each  class.  It  consists  of  an  image  that  conveys  the  computed 
probabilities  for  all  pixels.  This  may  be  a  grayscale  image  such  that  black  and  white 
represents  zero  and  unity  probabilities  respectively.  The  classified  image  consists  of  a 
number  of  uniform  patches  identified  by  a  particular  color  that  associates  each  patch  with 
a  class. 


b.  Spectral  A  ngle  Mapper 

The  spectral  angle  mapper  (SAM)  classifies  pixels  based  on  their  spectral 
similarity  to  reference  spectra  or  'end  member'  spectra.  A  computer  with  appropriate 
software  computes  spectral  similarity  by  treating  pixels  as  vectors  in  digital  number 
space.  This  space  is  consequently  band  intensity  space.  A  vector  can  be  created  by 
establishing  its  tail  at  the  origin  and  terminating  at  its  representative  pixel.  The  length  of 
the  vector's  projection  onto  a  particular  band's  axis  corresponds  to  the  pixel's  brightness 
in  that  band. 

The  SAM  determines  spectral  similarity  by  computing  the  angle  (spectral 
angle)  between  vectors  that  represent  the  reference  spectrum  and  the  unclassified  scene 
spectrum.  Figure  4.4  illustrates  the  idea  of  spectral  angle. 

When  the  reference  spectrum  and  the  unclassified  scene  spectrum  are 
composed  of  many  pixels,  Equation  4.3  is  used  to  compute  the  spectral  angle. 


[Equation  4.3] 
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a  =  Spectral  angle 

t.  =  Unclassified  scene  spectrum 

r.  =  Reference  spectra  vector  components 

N  =  Number  of  bands  of  data 

Take  note  that  a  is  independent  of  vector  length.  This  means  that  the  SAM  does  not 
consider  pixel  brightness.  This  allows  comparison  between  library  spectra  and  the 
spectra  of  target  pixels  even  with  incomplete  knowledge  of  the  illumination 
characteristics  of  the  target  pixel  spectra  (Kruse,  1993). 

In  the  same  likeness  as  the  maximum  likelihood  classifier,  with  the  SAM 
the  user  can  input  a  spectral  angle  threshold.  The  computed  values  of  a  can  be  directly 
translated  to  the  0-255  scale  for  display  purposes  (Mather,  1999).  The  output  of  the  SAM 
is  a  classified  image  and  a  number  of  decision  rules  that  equals  the  number  of 
endmembers.  The  decision  rules  may  be  grayscale  images,  such  that  black  and  white 
represents  large  and  small  angular  deviations  from  the  reference  spectrum  vector 
respectively.  The  classified  image  consists  of  a  number  of  uniform  patches  identified  by 
a  particular  color  that  associates  each  patch  with  a  class. 

D.  Accuracy  Assessment 

In  order  to  evaluate  the  access  the  accuracy  of  the  classification  process,  the 
software  tool  computes  a  confusion  matrix.  The  confusion  matrix  is  a  labeled  table 
illustrating  the  results  of  a  comparison  between  a  reference  image  and  a  classified  image. 
Comparisons  are  made  on  a  pixel-by-pixel  basis.  The  confusion  matrix  gives  two  kinds 
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of  errors  for  each  informational  class.  They  are  called  omission  and  commission  errors. 
Omission  errors  quantify  the  portion  of  pixels  in  a  particular  reference  class  that  the 
classifier  assigned  to  other  classes.  Commission  errors  occur  when  a  particular  class  in 
the  classified  image  is  incorrectly  assigned  pixels  from  other  classes. 

The  confusion  matrix  also  expresses  consumer  and  user  accuracy.  Producer 
accuracy  quantifies  the  portion  of  pixels  from  a  particular  reference  class  that  were 
classified  as  belonging  to  that  same  class.  Consumer  accuracy  is  the  ratio  of  correctly 
classified  pixels  assigned  to  a  particular  class  to  the  total  number  of  pixels  classified 
(correctly  and  incorrectly)  as  that  class. 

A  percent  correct  and  a  coefficient  of  agreement,  Kappa,  are  additional 
parameters  computed  in  the  confusion  matrix.  Kappa  measures  the  agreement,  beyond 
chance,  between  two  maps,  taking  into  account  all  pixels  of  an  image  (Wilkie  and  Finn, 
1996).  The  confusion  matrix  is  a  square  matrix  with  rows  corresponding  to  reference 
classes  and  columns  corresponding  to  classified  image  classes.  Percentage  correct  is  the 
percentage  ratio  of  the  number  of  pixels  correctly  classified  to  the  number  of  pixels 
examined.  Percentage  correct  may  be  given  for  each  class  and  for  the  overall 
classification.  Since  the  confusion  matrix  is  a  square  matrix,  number  of  pixels  correctly 
classified  is  the  sum  of  the  diagonal  elements.  Kappa  expresses  the  agreement  between 
two  images.  Equations  4.4,  4.5,  and  4.6  show  how  elements  of  the  confusion  matrix  are 
computed. 

c 

Total  percentage  Correct  =  [Equation  4.4] 

e.f  =  Confusion  matrix  elements 
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P  =  Total  number  of  pixels 


c  =  Number  of  elasses 

c 

Class  pereentage  eorreet  =  —  [Equation  4.5] 

7=1 

Z  A,  -  Z  [pi^dpiyd] 

Kappa  =  — - -  [Equation  4.6] 

i,7=l 

p.j  =  The  probability  of  the  ^  pixel  in  image  x  being  the  pixel 
in  image  y 

c 

^ =  The  pereentage  eorreet 
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V. 


LANDSAT  7  AND  IKONOS 


A,  OPERATIONAL  DIFFERENCES 

The  output  of  the  Ikonos  and  Landsat  7  satellite  systems  differ  in  spatial  and 

speetral  resolution.  Ikonos  provides  4-meter  4-band  multispeetral  and  1 -meter 
panehromatie  imagery.  Landsat  7  provides  30-meter  multispeetral  (bands  1-5  and  7),  15- 
meter  panehromatie  (band  8),  and  60-meter  thermal  (band  6)  imagery. 

B,  LANDSAT  7 

1.  Enhanced  Thematic  Mapper  Plus  (ETM+) 

Landsat  7  was  built  by  Lockheed  Martin  Missiles  and  Space  and  was  launched  on 
April  15,  1999.  It  is  the  most  recently  successfully  launched  satellite  in  the  Landsat 
series.  Landsat  7  uses  the  Enhanced  Thematic  Mapper  Plus  (ETM+)  sensor,  built  by 
Raytheon  Santa  Barbara  Remote  Sensing,  to  measure  radiation  reflected  or  emitted  from 
the  earth.  The  sensor  measures  radiation  in  eight  bands  and  produces  better  resolution  in 
the  thermal  infrared  than  the  Thematic  Mapper(TM)  sensor  onboard  Landsat  4  and  5. 
The  eight  bands  detect  radiation  in  the  visible,  near  infrared,  short-wave,  and  thermal 
infrared  regions  of  the  electromagnetic  spectrum. 

The  ETM+  collects  radiation  using  a  scanning  mirror  that  sweeps  across  the 
direction  of  travel.  A  telescope  focuses  the  radiation  into  compensation  mirrors,  which 
then  sends  it  to  the  focal  planes.  Compensation  mirrors  are  used  to  correct  for  scan  lines. 
Scan  lines  is  the  compound  effect  of  along-track  orbital  motion  and  cross-track  scanning 
which  leads  to  significant  overlap  and  underlap  in  ground  coverage  between  successive 
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scans  (SDUH).  The  energy  is  then  partitioned.  Some  of  the  partitioned  energy  is  sent  to 
the  silieon  deteetors  for  bands  1-4  and  8.  The  rest  of  the  energy  is  sent  to  a  cold  focal 
plane  (maintained  at  91  if  )  that  eontains  the  detectors  for  bands  7,  5,  and  6. 


2.  Orbit 

Landsat  7  operates  in  a  sun  synchronous  orbit  at  an  altitude  of  705  km .  This 
allows  it  to  pass  over  the  same  point  on  earth  at  essentially  the  same  local  time  every  16 
days  (SDUM).  However  the  sun's  elevation  angle  may  change  significantly  over  this 
time.  The  production  of  a  useful  image  requires  consideration  of  dependencies  of 
material  refieetances  and  atmospherie  effects  on  the  sun's  elevation  angle.  The  elevation 
angle  made  between  a  plane  tangent  to  the  earth's  surface  and  a  line  leading  to  the  sun. 
Landsat  7  is  positioned  perpendicular  to  this  tangent  plane.  Images  can  be  produeed  most 
efficiently  at  sun  elevation  angles  above  30°  and  a  low  angle  limit  has  been  set  at  15  ° . 

3.  Radiometric  Calibration  and  Relative  Spectral  Response 

Two  spherical  integrating  sourees  were  used  to  ealibrate  the  ETM-H  prior  to 
launeh  (SDUH).  These  deviees  consist  of  a  hollow  sphere  that  has  a  small  hole  and  a 
uniformly  eoated  inner  surfaee.  A  radiation  source  is  plaeed  inside  the  sphere.  Radiation 
exits  the  sphere  through  the  small  hole  in  a  random  orientation  and  provides  a  radiance 
speetrum  covering  the  entire  range  of  the  deteetor.  Following  SDUH,  calibration 
reduction  on  all  bands  except  band  6  was  performed  using  Equations  5.1  and  5.2. 
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[Equation  5.1] 


Lx{b,s) 


\RSR{b,X)L^{s,X)dX 

\RSR{b,X)dX 


=Band  Weighted  Speetral  Radiances  for  band  b  and 
sphere  level  5 

RSR{b,  A)  =  Relative  Spectral  Response  from  band  b  at  A 
L^{s,A)=  Measured  Spectral  Radiance  of  Sphere  level  s  at  A 

Q{d,b,s)  =  G{d,b)L^{b,s)  +  B{d,b)  [Equation  5.2] 

Q{b,  d,s)  =  Quantized  detector  responses 
G{d,b)=  Gains 
B{d,b)  =  Biases 

In  Equation  5.2,  Q{d,b,s)  is  regressed  against 5)  . 

Band  6  is  the  thermal  band.  In  calibration  it  is  matched  with  the  Thematic 
Mapper  Calibrator  (TMC),  which  uses  selectable  blackbody  sources.  Also  since  the 
ETM+  itself  is  a  producer  of  thermal  radiation,  following  SDUH  contributions  of  this 
thermal  radiation  to  calibration  is  modeled  using  Equation  5.3. 

Qsc  -  Qsh  id)  =  G{d){L^  ^^  -  )  [Equation  5 .3] 

d^A,esh  “  ^  ^  j  id^sh  ~  ) 

Qscid)  =  Quantized  Response  of  Band  6  at  detector  d  to  the  scene 
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Qsh(d)=  Quantized  Response  of  Band  6  at  deteetor  d  to  the  shutter 


G(d)  =  Gain  of  the  Deteetor 

sc  =Speotral  Radianee  of  the  seene 
^Aesh^  Seene-Equivalent  Speetral  Blackbody  Radianee  of  the 
Shutter 

=  Blaekbody  Radianee  of  the  Shutter 
Lj  =Blaekbody  Radianee  of  the  eomponent  of  the  ETM+ 
aj  =  Emissivity  adjusted  view  factor  for  the  component  of  the 
ETM+ 


Linder  operation  calibration  is  performed  periodically  with  three  onboard 
calibration  devices.  They  are  the  Internal  Calibrator  (IC),  the  Partial  Aperture  Solar 
Calibrator  (PASC),  and  the  Pull  Aperture  Solar  Calibrator  (PASC).  The  IC  is  useful  for 
calibrating  all  ETM+  bands,  while  the  PASC  and  PASC  are  mainly  useful  for  the 
reflective  bands  (SDUH).  The  IC  consists  of  two  tungsten  lamps  as  blackbody  sources. 
The  PASC  consists  of  a  white  panel  that  can  be  positioned  in  front  of  the  ETM+  aperture. 
The  white  panel  diffusely  reflects  radiation  into  the  ETM+  full  aperture.  With  known 
surface  reflectance,  solar  irradiance  and  geometry  conditions,  this  device  behaves  as  an 
independent,  full  aperture  calibrator  (SDUH).  The  PASC  allows  the  ETM+  to  image  the 
sun  while  viewing  a  'dark  earth'.  The  PASC  is  attached  to  the  ETM+  sun  shade 
obscuring  a  small  portion  of  the  aperture.  Pigures  5.1-8  illustrate  the  spectral  responses 
of  the  ETM-i-  eight  bands. 
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The  ETM+  has  8hzY  radiometric  resolution.  This  is  optimized  by  adjusting  the 
gain  state.  The  gain  is  adjusted  to  low  or  high  for  high  or  low  surface  brightness 
respectively.  Figure  5.9  illustrates  this. 

Figure  5.10  illustrates  the  relative  location  of  the  ETM+  onboard  the  Fandsat  7. 

C.  IKONOS 

Ikonos  was  built  by  Fockheed  Martin  Commercial  Space  Systems  for  Space 
Imaging  and  was  launched  in  September  1999.  For  two  years  it  had  the  highest 
commercially  available  spatial  resolution  of  1  meter  panchromatic  and  4  meter  4  band 
multispectral.  It  has  a  nominal  swath  width  of  Wkm  at  nadir  (sensing  vertically 
downward)  allowing  for  enhanced  revisit  rate  and  stereo  capabilities.  The  sensor  was 
built  by  Eastman  Kodak  Company. 

1.  Sensor  System 

The  sensor  system  consists  of  an  optical  telescope  assembly,  a  focal  plane  unit,  a 
digital  processing  unit,  and  a  power  supply  unit.  The  camera  telescope  has  the  equivalent 
resolving  power  of  a  10,000  mm  telephoto  lens  ("Figures  of  Merit").  It  consists  of  five 
curved  mirrors  and  two  fiat  mirrors  that  allow  for  a  reduced  telescope  length  while 
focusing  radiation  onto  the  focal  plane.  The  primary  mirror  has  a  honeycomb  glass  core 
covered  with  a  thin  highly  polished  mirror. 

The  focal  plane  unit  contains  the  sensor  arrays.  The  sensor  arrays  consist  of 
separate  CCD  arrays  for  simultaneously  collecting  panchromatic  and  multispectral 
imagery.  Filters  are  used  to  isolate  the  red,  green,  blue,  and  near  infrared  portion  of  the 
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collected  radiation.  The  digital  data  from  the  CCD  array  eleetronies  is  eventually 
compressed  in  the  digital  proeessing  unit.  The  digital  processing  unit  compresses  the 
digital  image  fdes  from  llhtoper  pixel  (bbp)  data  to  an  average  value  of  2.6 bbp  at  a 
speed  of  1 15  million  pixels  per  seeond  ("Figures  of  Merit"). 

The  Power  supply  system  uses  hybrid  DC-DC  eonverters  and  electronic  filters  to 
convert  high  voltage  direet  current  (DC)  to  a  relatively  low  voltage  DC.  The  low  voltage 
DC  is  used  to  power  the  camera's  sub-systems.  Total  eamera  power  consumption  is 
350  watts  .  Figure  5.11  in  Appendix  A  provides  additional  specifications  for  the  Ikonos 
system. 


2,  Radiometric  Calibration  and  Relative  Spectral  Response 

a.  Relative  Spectral  Response 

Ikonos  spectral  band  characteristics  are  given  in  Table  5.1  on  the 
following  page. 

Figures  5.12  and  5.13  illustrate  the  relative  spectral  responses  for 
panchromatic  and  4-band  multispectral  sensing. 

The  Ikonos  spectral  responses  were  determined  analytically  by  combining 
the  spectral  reflectivity  data  of  optical  witness  samples  and  the  measures  net  quantum 
efficiency  of  the  spectral  filter/detector  assemblies  (Space  Imaging).  Equation  5.4  is 
integrated  over  the  wavelength  interval  of  350-1035  nm  for  computing  Ikonos  in-band 
radiance. 
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L,=^L{X)R,^{X)dX 


[Equation  5.4] 


L^.  =  In-band  radiance  at  the  sensor  aperture  for  speetral  band  k 
L{A)  =  Spectral  radiance  at  the  sensor  aperture 

(A)  =  Peak-normalized  speetral  response  for  speetral  band  k 

Ikonos  imagery  is  ealibrated  to  at-aperture  in-band  radiance  by  resealing 
its  transmitted  data.  The  DN  values  are  converted  to  in-band  radiance  with  Equation  5.5 
and  the  ealibration  eoefficients  from  Table  5.2. 

E.  j  ^  =  DN-  ji^*  [CalCoefi^  ]”'  [Equation  5 .5] 

i,j,k  =  Indices  representing  image  pixel  i,j  in  speetral  band  k 
A  j  k~  In-band  radianee  at  the  sensor  aperture  for  pixel  i,  j  in 
speetral  band  k 

CalCoefi^  =  In-band  radianee  calibration  coefficient  for  spectral 
band  k 

E)A^.  ^  =  Image  product  digital  numbers  pixel  i,j  for  spectral 
band  k 

The  ealibration  coeffieients  for  eaeh  spectral  band  k  were  eomputed  using  Equation  5.6 
(Proe.  SPIE  Vol.  4381). 
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CalCoef,  =  j  L(X)R,  {l)L-'dX 


[Equation  5.6] 


=  Gain  Coefficient  for  spectral  band  k 

Aopt  =  Clear  aperture  area  of  the  Ikonos  telescope 
=  Ikonos  IFOV  in  spectral  band  k 
L^.  =  In-band  solar  radiance  at  the  sensor  aperture  for  spectral  band 
k 

L{X)  =  Solar  spectral  irradiance  at  the  sensor  aperture 

Rk{A)  =  Spectral  responsivity  for  spectral  band  k 
/I  =  Wavelength 
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VI.  MULTISPECTRAL  ANALYSIS 


A,  DESCRIPTION  OF  THE  STUDY  AREA 

This  thesis  examined  the  Elkhorn  Slough  and  its  surrounding  wetlands.  The 
Elkhorn  Slough  (ES)  sits  on  the  eentral  California  eoast  about  100  miles  south  of  San 
Eraneiseo  (Silberstein  and  Campbell,  1989).  Its  waterway  is  a  narrow  winding  body  that 
measures  700  ft  wide  and  25  /tdeep  at  its  mouth  in  the  town  of  Moss  Eanding.  The  ES 
is  supplied  freshwater  from  winter  rains  and  salt  water  from  the  sea.  It  is  composed  of  a 
number  of  habitats  including  the  waterway,  mudflats,  salt  marsh,  and  the  uplands,  which 
mainly  consists  of  taller  plants. 

Eocal  industry  includes  National  Refractories,  which  produce  heat-resistant  bricks 
by  extracting  magnesium  from  seawater,  and  Pacific  Gas  and  Electric,  which  produces 
electricity.  Conservation  of  the  ES  is  currently  being  undertaken  by  the  Elkhorn  Slough 
Eoundation  (ESE). 

B,  COLLECTING  DATA 

Ikonos  satellite  imagery  of  the  ES  its  surrounding  wetlands  was  collected  for 
October  23,  2000.  The  other  input  for  the  study  was  a  classification  map  provided  by  the 
ESE.  The  vector  map  defined  geo-referenced  polygons  that  define  qualitative  attributes 
of  the  primary  regions-of-interest  (ROIs).  The  vector  format  file  was  in  standard 
ARCVIEW  shape  file  format.  The  map  was  the  product  of  a  compilation  of  information 
from  aerial  photos  and  Global  Positioning  System  (GPS)  assisted  field  observations.  The 
ESP  used  a  50%  thumb  rule  to  qualitatively  define  each  ROT  The  thumb  rule  stated  that 
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each  region  was  defined  to  be  the  material  which  comprised  at  least  50%  of  it.  Tables  6. 1 
and  6.2  in  Appendix  B  provides  the  nomenclature  for  the  ROIs. 


C.  PROCESSING  THE  DATA 

The  bulk  of  the  analysis  time  was  spent  resolving  the  ROIs  derived  from  the 
above  noted  class  map.  A  fairly  standard  process  was  used  in  the  Environment  for 
Visualizing  Images  (ENVI).  A  Principle  Components  (PC)  transform  was  applied  to  the 
four  band  data,  which  facilitated  visualizing  the  data. 

The  primary  problem  with  the  initial  set  of  ROIs  was  that  most  reflected  a 
significant  mixture  of  regions  (e.g.  soil  and  vegetation),  and  these  needed  to  be  resolved 
or  purified. 

D,  TRAINING  DATA 

In  order  to  use  the  ROIs  as  training  data  for  the  classification  process,  it  was 
necessary  to  refine  the  ROIs  spectrally  to  ensure  that  each  ROTs  spectral  information 
uniquely  defined  it.  Analysis  of  the  true  and  false  color  images  and  using  ENVTs  N- 
Dimensional  Visualizer  (NDV)  provided  the  means  to  refine  the  spectral  data. 

The  NDV  produces  a  color-coded  plot  of  pixels  in  Digital  Number  (DN)  space 
and  allows  user  directed  or  automated  rotation  of  user  perspective  through  various 
dimensions.  The  NDV  displays  a  two  dimensional  slice  of  DN  space  at  a  time.  The  DN 
space  had  four  dimensions,  one  for  each  spectral  band  of  the  Ikonos  imagery,  or  three 
dimensions  in  PC  space.  The  NDV  also  affords  the  user  the  ability  to  color  code  pixels 
during  visualization  and  to  export  these  color-coded  pixels  back  and  forth  between  the 
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NDV  and  the  image.  This  is  useful  for  isolating  pixels  during  visualization  and  for 
providing  a  eonvenient  way  to  link  speetral  and  spatial  information. 

In  the  first  step  of  the  refinement  proeess,  all  ROIs  were  individually  analyzed 
with  the  NDV  using  the  raw  data  and  the  PCA  transformed  data.  This  was  awkward 
beeause  of  the  large  number  of  pixels,  and  a  laek  of  well  defined  elusters  in  DN  spaee. 
This  was  addressed  by  defining  small  ROIs  whieh  eould  be  used  as  referenee  elasses. 
The  reference  classes  were  obtained  by  defining  spectrally  pure  field  vegetation,  sand, 
soil,  trees,  concrete,  and  water.  Also  a  shade  test  region  was  created  for  calibration  of  the 
NDV,  and  for  providing  a  means  for  determining  the  direction  of  increasing  brightness 
once  rotation  in  DN  space  has  started.  Figure  6.1  illustrates  the  reference  spectra  in  DN 
space.  Although  the  axes  are  not  seen,  the  horizontal  axis  is  the  near  infrared  band  and 
the  vertical  axis  is  the  green  band.  A  third  band,  which  is  also  unseen,  is  perpendicular  to 
this  page;  it  is  the  red  band.  Each  ROI  and  the  test  regions  were  then  simultaneously 
visualized  with  the  NDV.  Figure  6.2  shows  the  test  regions  and  the  Eucalyptus  ROI  in 
the  NDV.  Figure  6.2  has  the  same  band  dimensions  as  Figure  6.1. 

Reference  classes  that  are  undesirable  to  include  in  the  training  data  for  a 
particular  ROI,  are  termed  excludable.  After  rotation  through  all  dimensions,  all  pixels 
collocated  with  and  in  sufficiently  close  proximity  to  excludable  reference  classes  were 
color-coded  white.  Color  coding  pixels  white  means  that  when  classes  are  exported  from 
the  NDV  to  the  image,  white  classes  are  excluded  from  exportation.  The  result  of  this 
process  was  a  new  set  of  secondary  ROIs  that  was  made  more  spectrally  pure  by 
eliminating  excludable  sub-classes. 
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In  the  second  step  of  the  refinement  process,  the  objective  was  to  ensure  that  the 
first  step  did  not  exclude  too  much  of  the  primary  ROI  data.  Excluding  too  much  of  the 
primary  ROI  data  could  produce  training  data  that  underestimated  the  variation  and 
overestimated  the  distinctiveness  of  the  region  it  was  supposed  to  represent.  The 
secondary  ROIs  were  again  individually  analyzed  with  the  NDV  using  the  raw  data  and 
the  PCA  transformed  data.  The  secondary  ROIs  were  repeatedly  exported  to  and  from 
the  NDV  and  a  false  or  true  color  image. 

E,  CLASSIFICATION 

1.  Spectral  Angle  Mapper 

The  SAM  was  applied  to  the  raw  data  using  various  angular  thresholds.  The 
SAM  used  both  the  refined  and  unrefined  ROIs  individually  as  training  data.  The  angular 
thresholds  used  were  0.075,  0.100,  0.150,  0.250,  and  0.300  radians.  A  classified  image 
and  a  set  of  rule  images  were  produced  for  each  angular  threshold. 

2.  Maximum  Likelihood 

The  maximum  likelihood  classifier  (MAXL)  was  applied  to  the  raw  data  using 
various  probability  thresholds.  The  MAXL  used  both  the  refined  and  unrefined  ROIs  as 
training  data.  The  probability  thresholds  were  90%  and  85%.  A  set  of  rule  images,  one 
for  each  region,  and  a  classified  image  was  produced  for  each  probability  threshold. 
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F.  CONFUSION  MATRIX 

A  confusion  matrix  was  computed  using  the  ENVI  result  for  all  elassified  images 
using  various  eombinations  of  the  refined  and  unrefined  training  data  as  ground  truth 
data.  Tables  6.2a-f  and  6.3a-f  provide  the  eonfusion  matriees  for  the  SAM  and  MAXL  at 
0.250  radians  and  85%  threshold  respeetively.  The  refined  data  was  used  as  training  and 
the  unrefined  data  was  used  as  ground  truth  data.  Table  6.1  in  the  next  seetion  should  be 
referred  to  for  illustration  of  the  various  eombinations. 


G.  RESULTS 

Table  6.1  illustrates  the  results  of  the  various  olassifieation  sehemes  and 


eonfusion  matrix  eomputation  sehemes. 


Classifieation 

Soheme 

Threshold 

Ground 
Truth  Data 

Training 

Data 

Overall 

Aoouraoy 

SAM 

0.075  Radians 

Refined 

Refined 

15.25% 

SAM 

0.100  Radians 

Refined 

Refined 

16.59% 

SAM 

0.150  Radians 

Refined 

Refined 

17.12% 

SAM 

0.250  Radians 

Refined 

Refined 

29.78% 

SAM 

0.250  Radians 

Unrefined 

Refined 

17.16% 

SAM 

0.250  Radians 

Unrefined 

Unrefined 

10.70% 

SAM 

0.300  Radians 

Refined 

Refined 

17.16% 

MAXE 

85% 

Refined 

Refined 

14.27% 

MAXE 

85% 

Unrefined 

Refined 

8.22% 

MAXE 

85% 

Unrefined 

Unrefined 

3.71% 

MAXE 

90% 

Refined 

Refined 

6.40% 

Table  6.1  Results  of  Classifieation  Sehemes 


The  overall  aeeuraeies  for  all  olassifieation  sehemes  were  low.  Generally  there  was  a 
marked  differenoe  between  the  results  produoed  by  the  two  olassifieation  sehemes. 
Overall  the  SAM  produoed  olassifieation  results  in  oloser  agreement  to  the  ESF  olass 
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map.  However,  from  visual  inspection,  the  rule  images  produced  by  the  MAXL  provide 
superior  utility  for  informational  reference.  For  most  of  the  rule  images,  the  MAXL 
resulted  in  far  more  suppression  of  regions  that  did  not  belong  to  the  rule's  class  than  did 
the  SAM.  The  classified  images  for  the  SAM  at  0.250  and  the  MAXL  at  85%,  using  the 
refined  data  as  training  data  and  ground  truth  data,  produced  the  best  results  and  are 
presented  in  Figure  6.44  and  Figure  6.43  respectively.  In  comparing  the  overall  accuracy 
of  the  SAM  at  0.250  radians  threshold,  using  the  refined  data  versus  the  unrefined  data 
for  both  training  and  ground  truth  data,  the  refining  process  was  responsible  for  a 
278.31%  magnification  in  accuracy.  Similarly,  the  refining  process  using  the  MAXL  at 
85%  probability  threshold  resulted  in  a  384.63%  magnification  in  accuracy. 

Figures  6.41  and  6.42  are  the  true  and  false  color  images  respectively.  Figures 
6.3-21  and  6.22-40  provide  the  rule  images  for  the  SAM  at  0.250  radians  threshold  and 
MAXL  at  85%  probability  threshold  respectively. 
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VII.  SUMMARY  AND  CONCLUSION 


The  Ikonos  satellite  was  used  to  collect  high  spatial  resolution  multispectral 
imagery  of  the  Elkhorn  Slough  was  collected  on  October  23,  2001.  Training  data  was 
obtained  from  the  Elkhorn  Slough  Eoundation.  A  rigorous  spectral  refining  process  was 
applied  to  the  raw  and  PCA  transformed  data.  The  effectiveness  of  using  the  training  data 
for  scene  classification  was  analyzed  by  applying  the  Spectral  Angle  Mapper  and  the 
Maximum  Eikelihood  classifiers.  Thresholds  were  varied,  and  the  classifiers  used  both 
the  refined  and  unrefined  training  data  separately.  Confusion  matrices  were  computed 
using  various  combinations  of  refined  and  unrefined  data  as  training  data  and/or  ground 
truth  data. 

Table  6.1  from  the  previous  chapter  demonstrated  the  effectiveness  of  using  high 
resolution  4-band  multispectral  data  to  perform  supervised  scene  classification  as  asserted 
in  Chapter  I.  Although  overall  accuracy  was  low  using  all  classification  schemes,  other 
results  were  promising.  The  results  indicated  that  up  to  a  385%  improvement  in  accuracy 
can  be  achieved  using  the  spectral  refining  process  described  in  Chapter  VI.  The  385% 
improvement  represents  the  amplification  of  accuracy  as  a  percentage.  The  Maximum 
Eikelihood  classification  result  produced  the  least  overall  accuracy. 

The  nature  of  the  data  provided  by  the  Elkhorn  Slough  Eoundation  was  in  part 
responsible  for  the  low  overall  accuracies.  Some  Regions-of-Interest  were  very  large  and 
too  spectrally  diverse  to  produce  high  classification  accuracy.  Eor  example,  reducing  the 
number  of  spectral  sub-classes  in  the  Tidal  Wetland  region  as  far  as  possible  without 
loosing  its  overall  spectral  identity  still  resulted  in  Tidal  Wetland  being  too  spectrally 
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diverse.  A  similar  effeet  oeeurred  with  some  other  regions.  Low  aceuraeies  also  eome  as 
a  result  of  the  thumb  rule  the  Elkhorn  Slough  Foundation  used  to  define  their  Regions-of- 
Interest.  This  thumb  rule  unavoidably  produced  a  large  number  of  heterogeneities  in 
most  Regions-of-Interest. 
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APPENDIX  A.  GENERAL  FIGURES  AND  TABLES 


Electromagnetic  Spectrum 


wavelength 


Gamma 

rays 


X  -  rays 


(1  nm)  10 


Ultraviolet 
*-  Visible: 


Infrared 

(IR) 


(1m)  10 


Microwaves 


frequerx^ 


/'  Blue  0.4  -  0.5  nm 
Green  0.5  -  0.6  nm 
Red  0.6  -  0.7  urn 

INFRARED  (IR) 
TERMINOLOGY 

Photographic  IR  0.7  -  0.9  pm 

Very  near  IR  0.7  - 1 .0  pm 

Reflected  IR  0.7  -  3.0  pm 

Near  IR  0.7  -  3.0  pm 

Thermal  IR  3.0  - 1 ,000  pm 
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Fig.  2.2  Atmospheric  Windows  (Lillesand  and  Kiefer,  1979) 
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Albedo  (%) 


A  (/<m) 


Fig.  2.3  Spectral  Albedos  (Rees,  1990) 
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Wavelength  nm 

0.4  0.6  0.8  1.0 


Fig.  2.4  Depth  Penetration  of  radiation  into  Clear  Water  (Gibson,  2000) 
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SPECTRAL  RADIANT  EMITTANCE 


Fig.  2.5  Spectral  Emittance  of  Three  Black  Body  Sources  (Lo,  1986) 
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Normalized 

intensity 


Fig.  2.6  Fraunhofer  Diffraction  Pattern  for  Circular  Aperture  (Karim,  1990) 
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Component 

Percentage 

N2 

78.08 

O2 

20.94 

Ar 

0.93 

CO2 

0.0314 

O3 

0.00000004 

Table  2.1  Gaseous  Components  of  the  Atmosphere  (Gibson,  2000) 
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WAVELENGTH,  MICROMETERS 

Fig.  2.7  Solar  Irradiance  With  Atmospheric  Gas  Attenuation  (Seyrafi  and 

Hovanessian,  1993) 
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Fig.  2.8  Aerosol  Number  Density  (Seyrafi  and  Hovanessian,  1993) 
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Fig.  2.9  Beam  Wander  of  a  ray  of  Radiation  eaused  by  Turbulenee  (Seyrafi  and 


Hovanessian,  1993) 
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Fig.  3.1  Operational  Altitudes  of  various  remote  sensing  systems  (Gibson,  2000) 
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Flight  lirt« 


Fig.  3.2  Strip  Camera  Image  Colleetion  (Lillesand  and  Kiefer,  1987) 
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Fig.  3.3  Panoramic  Camera  Image  Colleetion  (Lillesand  and  Kiefer,  1987) 
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Fig.  3.4  Pushbroom  Scanning  (Gibson,  2000) 
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Wavelength  In  Micrometers 
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MULTISPECTRAL 

WIDE  BANDWIDTH 

MODERATE  SPECTRAL  RESOLUTION 


A 


Spectral  identification  of  major  features, 
i.e.,  trees,  grass,  roads 


HYPERSPECTRAL 
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HIGH  SPECTRAL  RESOLUTION 


Spectral  discrimination  by  species, 
materials,  and  environmental  conditions 


ULTRASPECTRAL 

VERY  NARROW  BANDWIDTH 
VERY  HIGH  SPECTRAL  RESOLUTION 


Identification  and  discrimination  of  subtle 
spectral  details  of  materials,  vapors,  and 
aerosols 


Fig.  3.4  Relative  Differenee  in  Speetral  Resolution  among  Speetral  Scanning 


Systems  (Multispectral  User's  Guide,  1995) 
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0.7  urn 


Fig.  4.1  Three  sets  of  multispeetral  data  with  different  speetral  resolutions 

(Gibson,  2000) 
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Fig.  4.3a  Spectral  Data  Before  Establishing  Principal  Components  (Wilkie  and 

Finn,  1996) 
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Fig.  4.3b  Spectral  Data  After  Establishing  Principal  Components  (Wilkie  and 

Finn,  1996) 
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Relative  Spectral  Response, 
RSR  (unitless) 


Landsat  Band  1 

Thematic  Mapper  Spectral  Response  . 


Fig.  5.1  Spectral  Response  of  the  ETM+  for  band  1  (SDUH) 
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Relative  Spectral  Response, 
RSR  (unitless) 


1 
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Fig.  5.2  Spectral  Response  of  the  ETM+  for  band  2  (SDUH) 


Landsat  Band  2 

Thematic  Mapper  Spectral  Response 
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Relative  Spectral  Response, 
RSR  (unitless) 


Landsat  Band  3 

Thematic  Mapper  Spectral  Response 


Wavelength,  \  (nm) 

Fig.  5.3  Spectral  Response  of  the  ETM+  for  band  3  (SDUH) 
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Relative  Spectral  Response, 
RSR  (unitless) 


Landsat  Band  4 

Thematic  Mapper  Spectral  Response  . 


Wavelength,  \  (nm) 

Fig.  5.4  Spectral  Response  of  the  ETM+  for  band  4  (SDUH) 
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Relative  Spectral  Response, 
RSR  (unitless) 


Landsat  Band  5 

Thematic  Mapper  Spectral  Response 


Wavelength,  \  (nm) 

Fig.  5.5  Spectral  Response  of  the  ETM+  for  band  5  (SDUH) 
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Relative  Spectral  Response, 
RSR  (unitless) 


Landsat  Band  6 

Thematic  Mapper  Spectral  Response 


0.2  - 
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Fig.  5.6  Spectral  Response  of  the  ETM+  for  band  6  (SDUH) 
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Relative  Spectral  Response, 
RSR  (unitless) 


Landsat  Band  7 

Thematic  Mapper  Spectral  Response 


Wavelength,  \(nm) 


Fig.  5.7  Spectral  Response  of  the  ETM+  for  band  7  (SDUH) 
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Relative  Spectral  Response, 
RSR  (unitless) 


Fig.  5.9  Reflective  Band  Gain  Dynamic  Ranges  (SDUH) 
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y 


+x 

Nadir 

Fig.  5.10  Landsat  7  with  ETM+  (SDUH) 
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Launch  Information 

Vehicle:  Lockheed  Martin  Athena  II 

Rocket  Location:  Vandenberg  Air  Force  Base,  California 

Resolution 

•  Ground  resolution  of  each  hand: 

1  -meter  panchromatic 
(nominal  at  <26deg  off  nadir) 

4-meter  multispectral 

•  The  ground  processing  software  has  the  capability  to  rapidly  process 
and  mosaic  imagery  collected  at  various  imaging  angles  so  as  to  create 
seamless  image  products  with  a  consistent  pixel  GSD. 

Revisit  frequency 

2.9  days  at  1 -meter  resolution;  1.5  days  at  1.5-meter  resolution.  These 
values  are  for  targets  at  40  degrees  latitude. 

Spectral  Bands 

Panchromatic: 

0.45  -  0.90  microns 

Multispectral: 

#1:  Blue  0.45 -0.52 

#2:  Green  0.52  -  0.60 

#3:  Red  0.63  -0.69 

#4:  Near  IR  0.76 -0.90 

(same  as  Landsat  4&5  TM  Bands  #1-4) 

Swath  Widths  & 
Scene  Sizes 

Nominal  swath  width: 

•  1 1  km  at  nadir 

Areas  of  interest: 

•  a  single  image  at  1 1  km  x  1 1  km 

•  strips  of  1 1  km  x  100  km  up  to  11  km  x  1000  km 

•  image  mosaics  of  up  to  12,000  sq.  km 

•  up  to  two  10,000  square  kilometer  areas  in  a  single  pass  within  a 
region 

Metric  Accuracy 

•12-meter  horizontal  and  10-meter  vertical  accuracy  with  no  control 

•  2-meter  horizontal  and  3 -meter  vertical  accuracy  with  GCP  control 

•  These  are  specified  as  90%  CE  (circular  error)  for  the  horizontal  and 
90%  LE  (linear  error)  for  the  vertical 

Orbital  Information 

Altitude: 

681  kilometers  /  423  miles 

Inclination  angle:  98.1  deg 

Speed:  4  miles  per  second  /  7  kilometers  per  second 

Descending  nodal  crossing  time:  10-11  a.m. 

Revisit  frequency:  3  days  at  1  -meter  resolution;  1.5  days  at  1.5  meter 
resolution 

Orhit  time:  98  minutes 

Orhit  type:  sun-synchronous 

Viewing  angle:  Agile  spacecraft  -  in-track  and  cross-track  pointing 
Weight:  1600  pounds 

Fig.  5.11  Additional  Ikonos  specifications  (Space  Imaging  SE-REF-016  Rev.  A) 
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Fig.  5.13  Multispectral  Relative  Spectral  Response  (Space  Imaging) 
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Band 

Lower 

50% 

(nm) 

Upper 

50% 

(nm) 

Bandwidth 

(nm) 

Center 

(nm) 

Pan 

525.8 

928.5 

403 

727.1 

MS-1 

(Blue) 

444.7 

516.0 

71.3 

480.3 

MS-2 

(Green) 

506 .4 

595.0 

88.6 

550.7 

MS-3 

(Red) 

631 .9 

697.7 

65.8 

664.8 

MS^ 

(VNIR) 

757.3 

852.7 

95.4 

805.0 

Note:  1.  Spectral  Bandwidths  are  Full-Width  at  Flalf-Max 

Table  5.1  IKONOS  Spectral  Band  Characteristics  (Space  Imaging) 
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IKONOS  Radiometric  Calibratioi 

7  Coefficients  for  11  bit  products  [mW/(cm^*sr*DN)] 

Production  Date 

Blue 

Green 

Red 

NIR 

1 

pre  2/22/01 

633 

649 

840 

746 

post  2/22/01 

728 

727 

949 

843 

IKONOS  Radiometric  Calibratio 

n  Coefficients  for  8  bit  products  [mW/(cm^*sr*DN)] 

Production  Date 

Blue 

Green 

Red 

NIR 

pre  2/22/01 

79 

81 

105 

93 

post  2/22/01 

91 

91 

119 

105 

Table  5.2  Calibration  Coefficients  for  8  and  1 1  bit  IKONOS  data  (Space  Imaging) 
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Fig.  6.1  NDV  Plot  of  Reference  Classes 
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Fig.  6.2  NDV  Plot  of  Reference  Class  with  Eucalyptus  ROI  (Sienna) 
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APPENDIX  B,  FIGURES  AND  TABLES  OF  RESULTS 


Figure  6.4  SAM  Conifer  Rule  Image 
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Figure  6.5  SAM  Cultivated  Fields  Rule  Image 
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Figure  6.8  SAM  Eucalyptus  Rule  Image 
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Figure  6.10  SAM  Fresh  Water  Marsh  Rule  Image 
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i2 


Figure  6.13  SAM  Ponds  Rule  Image 
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Figure  6.16  SAM  Shadowed  Oak  Rule  Image 
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Sunlit  Oak  Rule  Image 


'8 


Figure  6.22  MAX.L.  Chaparral  Rule  Image 
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Figure  6.23  MAX.L.  Conifer  Rule  Image 
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Figure  6.24  MAX.L.  Cultivated  Fields  Rule  Image 
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Figure  6.25  MAX.L.  Development  Rule  Image 
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Figure  6.26  MAX.L.  Dry  Ponds  Rule  Image 
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Figure  6.27  MAX.L.  Eucalyptus  Rule  Image 
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Figure  6.28  MAX.L.  Fallow  Rule  Image 
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•>- 


Figure  6.29  MAX.L.  Fresh  Water  Marsh  Rule  Image 
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Figure  6.30  MAX.L.  Green  House  Rule  Image 
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Figure  6.31  MAX.L.  Maritime  Chaparral  Rule  Image 
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Figure  6.32  MAX.L.  Ponds  Rule  Image 
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Figure  6.33  MAX.L.  Riparian  Rule  Image 
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Figure  6.34  MAX.L.  Sage  Scrub  Rule  Image 
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Figure  6.36  MAX.L.  Slough  Water  Rule  Image 
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Figure  6.37  MAX.L.  Sunlit  Oak  Rule  Image 
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Figure  6.38  MAX.L.  Tidal  Wetlands  Rule  Image 
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Figure  6.39  MAX.L.  Turf  Rule  Image 
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Figure  6.40  MAX.L.  Dune  Serub  Rule  Image 
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Figure  6.41  True  Color  Image  of  the  Elkhorn  Slough 
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Figure  6.42  False  Color  Image  of  the  Flkhom  Slough 


132 


Figure  6.43  MAX.L.  (85%  Threshold)  Classified  Image 
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Figure  6.44  SAM  (.250  Radian  Threshold)  Classified  Image 
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Ground  Truth 

(Percent) 

Class 

CONIFER 

DEVELOPMENT 

EUCALYPTUS 

FALLOW 

FWM 

Unclassified 

0.00 

0.01 

0.00 

0.00 

0.00 

CONIFERS 

5.31 

0.46 

6.75 

0.10 

0.41 

DEVELOPMENT 

0.35 

3  6.77 

0.04 

1.80 

0.57 

EUCALYPTUS 

8.10 

0.32 

11.21 

0.09 

0.38 

FALLOW 

2.47 

8.12 

0.41 

42.99 

8.48 

FWH 

3.02 

3 . 15 

1.27 

3.02 

13 . 61 

GREENHOUSE 

0.48 

15.70 

0.09 

11.96 

0.71 

H.  CHAPARRAL 

2 . 14 

1.01 

1.46 

0.26 

2.27 

RIPARIAN 

1.75 

0.84 

1. 15 

0.29 

2 . 68 

SAGE  SCRUB 

12 . 13 

4.15 

9.56 

4.06 

8.52 

SLOUGH  WATER 

1.89 

1.96 

0.56 

0.07 

0.24 

TURF 

3.27 

1.26 

5. 68 

2.36 

7.62 

CULT.  FIELDS 

2 . 16 

2 . 16 

2.79 

1.83 

13 .46 

DUNE  SCRUB 

1.79 

1.85 

0.51 

4.75 

7.79 

C HAP P ARAL 

5.84 

1.77 

6.40 

0.70 

6.46 

SUNLIT  OAK 

15.65 

0.92 

22.42 

0.41 

3.92 

SHADOWED  OAK 

10.05 

0.99 

15.52 

0.21 

1.04 

DRY  PONDS 

5.47 

6.63 

0.96 

20. 62 

8.97 

TIDAL  WETLAND 

2.56 

2 . 15 

0.98 

2.42 

10.00 

PONDS 

15.59 

9.79 

12.25 

2.02 

2.86 

Total 

100.00 

100.00 

100.00 

100.00 

100.00 

Table  6.1a  SAM  Confusion  Matrix 


Ground  Truth 

(Percent) 

Class 

GREENHOUSE 

M.  CHAPARRAL 

RIPARIAN 

SAGE  SCRUB 

SLOUGH  WATER 

Unclassified 

0.00 

0.00 

0. 10 

0.00 

0.00 

CONIFERS 

0.02 

5.21 

2.29 

1.70 

0.00 

DEVELOPMENT 

28.24 

0.26 

0.90 

0.07 

2.39 

EUCALYPTUS 

0.03 

4.37 

1.97 

1.74 

0.00 

FALLOW 

13 .41 

1.31 

5.74 

4.06 

1.86 

FWM 

4.21 

3.51 

5.39 

3 . 17 

1.39 

GREENHOUSE 

2  6.57 

0.56 

0.85 

0.39 

1.06 

M.  CHAPARRAL 

0.09 

9.51 

5.64 

3.41 

0.05 

RIPARIAN 

0.09 

7.32 

4.55 

2.51 

0.11 

SAGE  SCRUB 

1.34 

13 .82 

13.87 

37.20 

1.75 

SLOUGH  WATER 

0.10 

0.29 

0.53 

0.63 

78. 17 

TURF 

0.45 

0.97 

5.34 

1.00 

0.88 

CULT.  FIELDS 

0.67 

2.44 

6.84 

2.21 

1.40 

DUNE  SCRUB 

2.03 

3 . 17 

4.21 

8.30 

0.84 

CHAPPARAL 

0.38 

13.73 

9.08 

4.71 

0.33 

SUNLIT  OAK 

0.35 

8.15 

7.83 

4.45 

0.09 

SHADOWED  OAK 

0.06 

9.74 

5. 15 

4.08 

0.02 

DRY  PONDS 

14.63 

3.48 

5.76 

4.38 

0.63 

TIDAL  WETLAND 

1.61 

6.82 

7. 11 

8.23 

0.80 

PONDS 

5.72 

5.35 

6.85 

7.76 

8.22 

Total 

100.00 

100.00 

100.00 

100.00 

100.00 

Table  6.1b  SAM  Confusion  Matrix 
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Class 

TURF 

Ground  Truth 

CULT.  FIELDS 

(Percent) 

DUNE  SCRUB 

CHAPPARAL 

SUNLIT  OAK 

Unclassified 

O.GG 

G.Gl 

0.00 

0.00 

0.00 

CONIFERS 

G.7G 

G.G4 

0.00 

2  .  16 

4.29 

DEVELOPMENT 

2 . 66 

4.53 

31.29 

0.00 

0.00 

EUCALYPTUS 

G.  65 

0.04 

0.00 

2.31 

9.98 

FALLOU 

2.24 

20.70 

26.43 

0.00 

0.00 

FWH 

8.29 

7.77 

1.21 

0.00 

0.00 

GREEN  HOUSE 

4.38 

5.22 

26.27 

0.00 

0.00 

H.  CHAPARRAL 

G.  67 

0.28 

0.08 

13 . 11 

0.05 

RIPARIAN 

G.  63 

0.84 

0.06 

16.43 

0.00 

SAGE  SCRUB 

3.87 

1.59 

4.65 

0.72 

1.45 

SLOUGHWATER 

G.  17 

0.08 

0.76 

0.00 

0.00 

TURF 

41.78 

11.62 

0.05 

0.00 

16.70 

CULT.  FIELDS 

11.34 

13 . 19 

0.27 

1.01 

0.00 

DUNE  SCRUB 

1.23 

1.41 

2.71 

0.00 

0.00 

CHAPPARAL 

1.75 

2.87 

0.23 

49.42 

3 . 15 

SUNLIT  OAK 

2.43 

2.49 

0.04 

11.38 

55.95 

SHADOWED  OAK 

1.78 

0.09 

0.00 

1.87 

7.96 

DRY  PONDS 

1G.76 

16.07 

3 . 14 

0.00 

0.00 

TIDAL  WETLAND 

1.87 

1.84 

1.40 

1 .59 

0.00 

PONDS 

2.77 

9.34 

1.41 

0.00 

0.47 

Total 

IGG.GG 

100.00 

100.00 

100.00 

100.00 

Table  6.1e  SAM  Confusion  Matrix 

Class 

SHADOWED  OAK 

Ground  Truth 

DRY  PONDS 

(Percent) 

TIDAL  WETLAND 

PONDS 

Total 

Unclassified 

0.00 

0.00 

0.00 

0.00 

0.01 

CONIFERS 

10.59 

0.20 

0.04 

0.35 

0.97 

DEVELOPMENT 

0.00 

2 . 61 

3 . 66 

4.32 

5.18 

EUCALYPTUS 

14.87 

0.23 

0.03 

0.22 

1.10 

FALLOW 

0.00 

28.80 

5.49 

8.09 

14.69 

FWM 

0.00 

5.23 

5.90 

3.86 

5.95 

GREENHOUSE 

0.00 

3.54 

4.54 

4.18 

5.18 

H.  CHAPARRAL 

0.00 

0.81 

0.28 

3 . 60 

1.47 

RIPARIAN 

0.00 

0.90 

0.42 

1.69 

1.48 

SAGE  SCRUB 

12 .37 

7.98 

15.41 

16. 63 

6.42 

SLOUGH  WATER 

1.20 

0. 12 

22.00 

7.80 

5.96 

TURF 

0.23 

0.41 

4.34 

2.57 

7.58 

CULT.  FIELDS 

0.00 

1.34 

5.61 

2.45 

8.42 

DUNE  SCRUB 

0.00 

4.56 

6.22 

4.52 

2.72 

CHAPPARAL 

1.09 

1.42 

1.81 

1.74 

3.83 

SUNLIT  OAK 

11.08 

0.26 

0.43 

0.64 

3.58 

SHADOWED  OAK 

32.01 

0.38 

0.11 

0.95 

2.04 

DRY  PONDS 

0.00 

31.45 

3.32 

4. 10 

10.87 

TIDAL  WETLAND 

0.00 

4.65 

6.50 

5.20 

3.33 

PONDS 

16.57 

5.11 

13.89 

27. 13 

9.22 

Total 

100.00 

100.00 

100.00 

100.00 

100. 0i 

Table  6. Id  SAM  Confusion  Matrix 
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Class 

Commission 

Omission 

Commission 

Omiss ion 

(Percent) 

(Percent) 

(Pixels) 

(Pixels) 

CONIFERS 

99.06 

94.69 

46212/46652 

7848/8288 

DEVELOPMENT 

64.65 

63 .23 

160694/248577 

151106/238989 

EUCALYPTUS 

52.36 

88.79 

27537/52593 

198494/223550 

FALLOW 

82.57 

57.01 

582014/704887 

162946/285819 

FWH 

96.81 

86.39 

276539/285664 

57899/67024 

GREENHOUSE 

96.64 

73.43 

240064/248417 

23083/3 1436 

H.  CHAPARRAL 

42.06 

90.49 

29576/70315 

387694/428433 

RIPARIAN 

89.45 

95.45 

63673/71182 

157497/165006 

SAGE  SCRUB 

94. 62 

62 .80 

291602/308195 

28008/44601 

SLOUGH  WATER 

59.40 

21.83 

169843/285916 

32420/148493 

TURF 

95.06 

58.22 

346007/363978 

25039/43010 

CULT.  FIELDS 

23.87 

86.81 

96521/404301 

2026209/2333989 

DUNE  SCRUB 

99.65 

97.29 

130167/130626 

16470/16929 

CHAPPARAL 

99.81 

50.58 

183407/183750 

351/694 

SUNLIT  OAK 

99.37 

44.05 

170869/171951 

852/1934 

SHADOWED  OAK 

97.57 

67.99 

95613/97992 

5052/7431 

DRY  PONDS 

99.79 

68.55 

520753/521836 

2361/3444 

TIDAL  WETLAND 

71.29 

93.50 

113856/159712 

659853/705709 

PONDS 

97.27 

72.87 

430305/442390 

32456/44541 

Table  6.1®  SAM  Confusion  Matrix 


Class 

Prod.  Acc. 

User  Acc. 

Prod.  Acc. 

User  Acc. 

(Percent) 

(Percent) 

(Pixels) 

(Pixels) 

CONIFERS 

5.31 

0.94 

440/8288 

440/46652 

DEVELOPMENT 

3  6.77 

35.35 

87883/238989 

87883/248577 

EUCALYPTUS 

11.21 

47.64 

25056/223550 

25056/52593 

FALLOW 

42.99 

17.43 

122873/285819 

122873/704887 

FWM 

13 . 61 

3 . 19 

9125/67024 

9125/285664 

GREENHOUSE 

2  6.57 

3.36 

8353/31436 

8353/248417 

M.  CHAPARRAL 

9.51 

57.94 

40739/428433 

40739/70315 

RIPARIAN 

4.55 

10.55 

7509/165006 

7509/71182 

SAGE  SCRUB 

37.20 

5.38 

16593/44601 

16593/308195 

SLOUGH  WATER 

78.17 

40.60 

116073/148493 

116073/285916 

TURF 

41.78 

4.94 

17971/43010 

17971/363978 

CULT.  FIELDS 

13 . 19 

76.13 

307780/2333989 

307780/404301 

DUNE  SCRUB 

2.71 

0.35 

459/16929 

459/130626 

CHAPPARAL 

49.42 

0.19 

343/694 

343/183750 

SUNLIT  OAK 

55.95 

0.63 

1082/1934 

1082/171951 

SHADOWED  OAK 

32.01 

2.43 

2379/7431 

2379/97992 

DRY  PONDS 

3  1.4 

0.21 

1083/3444 

1083/521836 

TIDAL  WETLAND 

6.50 

28.71 

45856/705709 

45856/159712 

PONDS 

2  7. 13 

2.73 

12085/44541 

12085/442390 

Table  6.  If  SAM  Confusion  Matrix 
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Ground  Truth 

(Percent) 

Class 

CONIFER 

DEVELOPMENT 

EUCALYPTUS 

FALLOW 

FWM 

Unclassified 

96 . 14 

83 . 18 

98.01 

80.29 

96.09 

CONIFERS 

0. 00 

0.00 

0.00 

0.00 

0.00 

DEVELOPMENT 

0.01 

7.62 

0.01 

0.09 

0.07 

EUCALYPTUS 

0.00 

0.00 

0.00 

0.00 

0.00 

FALLOW 

0.49 

1.74 

0.03 

16.26 

1.60 

FWH 

0.00 

0.00 

0.00 

0.00 

0.00 

GREENHOUSE 

0.29 

3.31 

0.04 

0.20 

0. 18 

H.  CHAPARRAL 

0.00 

0.00 

0.00 

0.00 

0.00 

RIPARIAN 

0.00 

0.00 

0.00 

0.00 

0.00 

SAGE  SCRUB 

0.00 

0.00 

0.00 

0.00 

0.00 

SLOUGH  WATER 

0.57 

0.01 

0.01 

0.00 

0.00 

TURF 

0.53 

0.36 

0.03 

0.2  6 

0.67 

CULT.  FIELDS 

0. 11 

0.50 

0.05 

0.3  6 

0.50 

DUNE  SCRUB 

0.06 

0.25 

0.00 

0.08 

0. 17 

C HAP P ARAL 

0.08 

0.01 

0.24 

0.00 

0.09 

SUNLIT  OAK 

0.34 

0.01 

0.51 

0.01 

0.03 

SHADOWED  OAK 

0.27 

0.00 

0.73 

0.00 

0.00 

PONDS 

0.98 

1.12 

0.08 

0.05 

0.36 

TIDAL  WETLAND 

0.02 

0.08 

0.04 

0.03 

0.03 

DRY  PONDS 

0. 11 

1.80 

0.20 

2.35 

0.19 

Total 

100.00 

100.00 

100.00 

100.00 

100.00 

Table  6.2a  MAXL  Confusion  Matrix 

Class 

GREENHOUSE 

Ground  Truth 

M.  CHAPARRAL 

(Percent) 

RIPARIAN 

SAGE  SCRUB 

SLOUGH  WATER 

Unclassified 

53 .27 

96.68 

96.73 

98.33 

21.95 

CONIFERS 

0.00 

0.00 

0.00 

0.00 

0.00 

DEVELOPMENT 

3 . 41 

0.04 

0.17 

0.01 

0.65 

EUCALYPTUS 

0.00 

0.00 

0.00 

0.00 

0.00 

FALLOW 

0.45 

0.06 

0.67 

0.07 

0.31 

FWM 

0.00 

0.00 

0.00 

0.00 

0.00 

GREENHOUSE 

41.28 

0.29 

0.29 

0.15 

0.38 

H.  CHAPARRAL 

0.00 

0.00 

0.00 

0.00 

0.00 

RIPARIAN 

0.00 

0.00 

0.00 

0.00 

0.00 

SAGE  SCRUB 

0.00 

0.01 

0.01 

0.10 

0.00 

SLOUGH  WATER 

0.01 

0.00 

0.00 

0.00 

71.22 

TURF 

0.  12 

0.05 

0.27 

0.02 

0.02 

CULT.  FIELDS 

0.51 

0.02 

0.56 

0.02 

0.09 

DUNE  SCRUB 

0.  17 

0.03 

0.07 

0.02 

0.09 

C HAP P ARAL 

0.00 

1.43 

0.05 

0.11 

0.00 

SUNLIT  OAK 

0.02 

0. 18 

0.06 

0.20 

0.00 

SHADOWED  OAK 

0.00 

0.19 

0.04 

0.13 

0.00 

PONDS 

0. 16 

0.06 

0.30 

0.06 

5.06 

TIDAL  WETLAND 

0.03 

0.01 

0.03 

0.01 

0.22 

DRY  PONDS 

0.58 

0.96 

0.77 

0.79 

0.00 

Total 

100.00 

100.00 

100.00 

100.00 

100.00 

Table  6.2b  MAXL  Confusion  Matrix 
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Ground  Truth 

(Percent) 

Class 

TURF 

CULT.  FIELDS 

DUNE  SCRUB 

CHAPPARAL 

SUNLIT  OAK 

Unclassified 

51. 18 

80.2  6 

92 . 10 

87. 18 

90.85 

CONIFERS 

0.00 

0.00 

0.00 

0.00 

0.00 

DEVELOP KENT 

0. 13 

0.38 

1.21 

0.00 

0.00 

EUCALYPTUS 

0.00 

0.00 

0.00 

0.00 

0.00 

FALLOW 

0.43 

3.04 

1.73 

0.00 

0.00 

FWH 

0.00 

0.00 

0.00 

0.00 

0.00 

GREENHOUSE 

3.55 

0.24 

3.41 

0.00 

0.00 

H.  CHAPARRAL 

0.00 

0.00 

0.00 

0.00 

0.00 

RIPARIAN 

0.00 

0.00 

0.00 

0.00 

0.00 

SAGE  SCRUB 

0.00 

0.00 

0.00 

0.00 

0.00 

SLOUGH  WATER 

0.00 

0.00 

0.43 

0.00 

0.00 

TURF 

42 . 19 

1.13 

0.00 

0.00 

0.00 

CULT.  FIELDS 

0.30 

7.40 

0.01 

0.00 

0.00 

DUNE  SCRUB 

0.05 

0.11 

0.90 

0.00 

0.00 

C HAP P ARAL 

0.00 

0.00 

0.00 

12.82 

0.00 

SUNLIT  OAK 

0.03 

0.03 

0.00 

0.00 

9.15 

SHADOWED  OAK 

0.00 

0.00 

0.00 

0.00 

0.00 

PONDS 

0. 17 

0.05 

0.02 

0.00 

0.00 

TIDAL  WETLAND 

0.03 

0.01 

0.19 

0.00 

0.00 

DRY  PONDS 

1.94 

7.35 

0.00 

0.00 

0.00 

Total 

100.00 

100.00 

100.00 

100.00 

100.00 

Table  6.2c  MAXL  Confusion  Matrix 

Ground  Truth 

(Percent) 

Class 

SHADOWED  OAK 

PONDS 

TIDAL  WETLAND 

DRY  PONDS 

Total 

Unclassified 

96.54 

73.93 

72.91 

81.  62 

80.37 

CONIFERS 

0.00 

0.00 

0.00 

0.00 

0.00 

DEVELOPMENT 

0.00 

1.02 

0.32 

0.52 

0.68 

EUCALYPTUS 

0.00 

0.00 

0.00 

0.00 

0.00 

FALLOW 

0.00 

0.83 

1.67 

4.82 

2.87 

FWH 

0.00 

0.00 

0.00 

0.00 

0.00 

GREENHOUSE 

0.00 

1. 10 

1.86 

0.32 

0.95 

H.  CHAPARRAL 

0.00 

0.00 

0.00 

0.00 

0.00 

RIPARIAN 

0.00 

0.00 

0.00 

0.00 

0.00 

SAGE  SCRUB 

0.00 

0.01 

0.01 

0.03 

0.00 

SLOUGH  WATER 

0.00 

0.30 

13 .97 

0.00 

4.26 

TURF 

0.00 

0.95 

0.07 

0.20 

1.01 

CULT.  FIELDS 

0.00 

0.20 

0. 10 

0.06 

3.70 

DUNE  SCRUB 

0.00 

0.48 

0.13 

0.03 

0.11 

CHAPPARAL 

0.00 

0.00 

0.00 

0.00 

0.15 

SUNLIT  OAK 

0.26 

0.02 

0.01 

0.00 

0.06 

SHADOWED  OAK 

3 .20 

0.01 

0.00 

0.00 

0.06 

PONDS 

0.00 

20.10 

7.24 

1.57 

1.52 

TIDAL  WETLAND 

0.00 

0.65 

1.44 

0.  12 

0.24 

DRY  PONDS 

0.00 

0.41 

0.28 

10.71 

4.01 

Total 

100.00 

100.00 

100.00 

100.00 

100.00 

Table  6. 2d  MAXL  Confusion  Matrix 
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Class 

Commission 

Omission 

Commission 

Omission 

(Percent) 

(Percent) 

(Pixels) 

(Pixels) 

CONIFERS 

0.00 

100.00 

0/0 

8288/8288 

DEVELOPMENT 

44. 61 

92.38 

14665/32871 

220783/238989 

EUCALYPTUS 

0.00 

100.00 

0/0 

223550/223550 

FALLOW 

66.  19 

83.74 

91018/137505 

239332/285819 

FWH 

100.00 

100.00 

2/2 

67024/67024 

GREENHOUSE 

71.46 

58.72 

32485/45462 

18459/3143  6 

H.  CHAPARRAL 

0.00 

100.00 

0/0 

428433/428433 

RIPARIAN 

0.00 

100.00 

0/0 

165006/165006 

SAGE  SCRUB 

74.72 

99.90 

133/178 

44556/44601 

SLOUGH  WATER 

48.32 

28.78 

98889/204643 

42739/148493 

TURF 

62 . 39 

57.81 

30095/48239 

24866/43010 

CULT.  FIELDS 

2.79 

92.60 

4960/177684 

2161265/2333989 

DUNE  SCRUB 

97.14 

99.10 

5171/5323 

16777/16929 

C HAP P ARAL 

98.73 

87.18 

6926/7015 

605/694 

SUNLIT  OAK 

94.29 

90.85 

2925/3102 

1757/1934 

SHADOWED  OAK 

91.75 

96.80 

2648/2886 

7193/7431 

PONDS 

87.73 

79.90 

64026/72977 

35590/44541 

TIDAL  WETLAND 

12  .  15 

98.56 

1405/11568 

695546/705709 

DRY  PONDS 

99.81 

89.29 

192095/192464 

3075/3444 

Table  6.2" 

MAXL  Confusion  Matrix 

Class 

Prod.  Acc. 

User  Acc. 

Prod.  Acc. 

User  Acc. 

(Percent) 

(Percent) 

(Pixels) 

(Pixels) 

CONIFERS 

0.00 

0.00 

0/8288 

0/0 

DEVELOPMENT 

7.62 

55.39 

18206/238989 

18206/32871 

EUCALYPTUS 

0.00 

0.00 

0/223550 

0/0 

FALLOW 

16.26 

33.81 

46487/285819 

46487/137505 

FWM 

0.00 

0.00 

0/67024 

0/2 

GREENHOUSE 

41.28 

28.54 

12977/31436 

12977/45462 

M.  CHAPARRAL 

0.00 

0.00 

0/428433 

0/0 

RIPARIAN 

0.00 

0.00 

0/165006 

0/0 

SAGE  SCRUB 

0.10 

25.28 

45/44601 

45/178 

SLOUGH  WATER 

71.22 

51.68 

105754/148493 

105754/204643 

TURF 

42 . 19 

37.61 

18144/43010 

18144/48239 

CULT.  FIELDS 

7.40 

97.21 

172724/2333989 

172724/177684 

DUNE  SCRUB 

0.90 

2.86 

152/16929 

152/5323 

CHAPPARAL 

12 . 82 

1.27 

89/694 

89/7015 

SUNLIT  OAK 

9.15 

5.71 

177/1934 

177/3102 

SHADOWED  OAK 

3.20 

8.25 

238/7431 

238/2886 

PONDS 

20. 10 

12 .27 

8951/44541 

8951/72977 

TIDAL  WETLAND 

1.44 

87.85 

10163/705709 

10163/11568 

DRY  PONDS 

10.71 

0.19 

369/3444 

369/192464 

Table  6.2f  MAXL  Confusion  Matrix 
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